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cesses  in  the  upper  atmosphere.  A  physical  model  that  includes  ex¬ 
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that  lead  to  emission  of  infrared  radiation  has  been  implemented  in 
a  computer  program  that  computes  infrared  radiances  for  an  earth’s 
limb  viewing  geometry.  The  nominal  spectral  region  of  this  study 
lies  between  five  and  twenty-five  micrometers  and  emphasis  is  plactj 
on  radiation  originating  at  altitudes  between  70  and  500  km.  An 
earlier  model  for  the  transport  of  infrared  radiation  in  molecular 
bands  with  Doppler  line  shape,  accurate  only  for  linear  molecules, 
has  been  extended  to  bands  of  polyatomic  molecules.  A  chemistry 
program  which  includes  the  effects  of  vertical  transport  by  eddy 
mixing  and  molecular  diffusion  has  been  developed  to  make  possible 
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SECTION  I 


INTRODUCTION 

The  research  reported  here  has  had  as  its  major  objective  the 
development  of  a  computer  program  to  simulate  the  natural  infrared 
radiance  background  of  the  earth's  upper  atmosphere.  The  nominal 
spectral  region  under  study  lies  between  five  and  twenty-five  micrometers 
and  emphasis  is  placed  on  radiation  originating  at  altitudes  between  70  and 
500  kilometers.  The  general  probrlem  area  is  of  interest  for  systems 
design,  military  surveillance  and  the  advancement  of  knowledge  about 
physical  processes  in  the  upper  atmosphere.  The  immediate  application 
of  this  work  will  be  to  aid  in  developing  optimum  infrared  background 
measurements  programs  and  in  interpreting  the  results  of  such  measure¬ 
ments. 

This  work  is  an  extension  of  the  study  of  Corbin,  et  al.  (1969)  who 
investigated  the  natural  infrared  background  of  the  earth  in  the  5  to  25  micrometer 
spectral  region,  with  the  goal  of  estimating  earth  limb  viewing  radiances  for 
tangent  heights  from  the  surface  to  500  km  altitude.  For  convenience,  their 
study  divided  the  atmosphere  into  two  regions  with  a  division  at  70  km.  Below 
70  km  the  atmosphere  was  assumed  to  be  in  thermal  equilibrium.  Above  70  km 
explicit  calculations  were  made  of  processes  which  excite  and  de-excite 
molecular  vibrational  and  rotational  levels  which  are  the  source  of  infrared 
radiation.  Their  study  concentrated  on  radiation  from  water  vapor,  carbon 
dioxide,  ozone,  nitric  oxide  and  nitrous  oxide,  which  are  the  principal  radiating 
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species  in  the  spectral  region  considered.  In  addition,  nitric  acid -was 
included  in  the  lower  atmosphere  work  and  estimates  were  made  of  radiation 
to  be  expected  from  particulate  matter  suspended  in  the  atmosphere. 

Corbin,  et  al.  (1969)  presented  models  for  the  lower  atmosphere  for 
a  wide  range  of  seasonal  and  latitudinal  conditions.  This  was  not  possible 
for  the  abundances  of  most  minor  neutral  species.  More  data  has  since  become 
available,  particularly  for  nitric  oxide  and  the  hydroxyl  radical,  but  at  present 
it  appears  that  the  best  means  of  estimating  abundances  of  important  infrared 
emitting  species  is  chemical  rate  equation  integrations  including  molecular 
diffusion  and  eddy  mixing.  Even  calculations  involving  transport  properties 
cannot  always  be  accepted  because  the  values  of  eddy  mixing  coefficients  are 
to  a  large  extent  only  informed  guesses  and  published  calculations  often  use 
outdated  rate  coefficients.  For  these  reasons,  one  of  the  tasks  of  the  current 
contract  has  been  to  develop  a  computer  program  with  which  to  determine 
diurnal  variations  in  abundances  of  minor  species,  to  investigate  the  effects 
of  changes  in  assumed  eddy  mixing  coefficients,  and  to  quickly  incorporate 
new  determinations  of  chemical  rate  coefficients. 

A  second  area  of  study  involved  in  improving  the  radiance  model 
includes  the  physical  processes  that  control  the  population  of  infrared  emitting 
states  of  atmospheric  molecules.  Except  for  the  pure  rotational  radiation  from 
molecules  such  as  water,  the  degree  of  excitation  of  vibrational  levels 
determines  the  radiation  from  infrared  emitting  molecules.  The  most  important 
mechanisms  are  collisional  excitation  and  de-excitation  and  absorption  and 
re-emission  of  electromagnetic  radiation. 

In  the  troposphere  and  lower  stratosphere,  collisional  processes  are 
rapid  enough  to  control  the  population  of  vibrational  levels.  Above  30  to  50  km, 
however,  collisional  excitation  becomes  less  efficient  and  radiative  processes 
become  important.  The  combined  effects  of  collisional  and  radiative  processes 
must  therefore  be  considered.  Below  an  altitude  of  about  90  km  molecular 
nitrogen  and  oxygen  are  the  most  important  collision  partners.  Above  that 


altitude,  atomic  oxygen  becomes  important,  both  in  exciting  nitric  oxide 
and  in  determining  molecular  oxygen  vibrational  temperatures.  Corbin, 
et  al.  (1969)  included  collisional  excitation  by  molecular  nitrogen  and  oxygen 
in  their  model  and  the  work  reported  here  adds  the  effects  of  atomic  oxygen 
as  a  collision  partner. 

In  determining  the  effects  of  radiation  on  the  populations  of  vibrational 
levels  of  infrared  active  molecules,  it  is  necessary  to  separate  the  radiation 
of  a  single  change  in  vibrational  quantum  numbers  from  the  rest  of  the 
radiation  field.  Corbin  et  al.  (1969)  did  this  by  assuming  a  Doppler  line  shape 
for  the  individual  rotational  lines  of  a  band  and  were  able  to  obtain  adequate 
numerical  approximations  for  radiative  transfer  functions  appropriate  to  single 
bands  of  linear  molecules.  This  approach  has  been  extended  to  water  vapor 
and  ozone  bands  lying  in  the  spectral  region  of  interest. 

The  following  sections  of  this  report  summarize  the  research  performed 
during  this  contract.  The  computational  algorithm  used  in  the  chemistry  program, 
relevant  chemical  reaction  rate  coefficients  and  currently  known  chemical 
abundances  in  the  upper  atmosphere  are  discussed  in  Section  II.  Relevant1 
collisional  and  radiative  excitation  mechanisms  are  reviewed  in  Section  HI, 
where  the  computational  model  for  predicting  high  altitude  infrared  radiance 
is  developed.  The  chemistry  and  infrared  radiance  programs  are  described 
and  listed  in  Section  IV.  Sample  results  for  White  Sands,  New  Mexico  are 
presented  in  Section  V,  where  limitations  of  the  model  and  suggestions  for 
future  v/ork  are  discussed. 


SECTION  II 


CHEMISTRY 


There  hu  **■  en  relatively  le w  measurements  of  the  variation  with  : 

»  ' 

altitude  of  infr«.  e«  emitting  species  at  altitudes  above  30  km,  the  limit  of 

i 

direct  sampling  by  balloons.  While  ozone  measurements  by  different 
workers  og  ee  reasonably  well  with  each  other  and  with  theoretical 
computations  at  altitudes  up  to  100  km,  and  nitric  oxide  profiles  are 
reasonably  well  known  between  about  80  and  150  km,  corresponding  data 
are  lacking  for  other  species  which  are  important.  There  is  in  some  cases 
indirect  data  available  which  when  combined  with  computational  models  can 

yield  limits  to  concentrationsi  Thus,  Anderson  (1971a)  has  been  able  to  put 

•  : 

an  upper  limit  cn  total  hydroxyl  radical  abundance  above  60  km  from  measure¬ 
ment  of  resonance  fluorescer.:e  and  with  this  result  computes  possible  water 

:  ;  I 

vapor  profiles  between  50  and  100  km  using  a  photochemical  model  that 
includes  vertical  transport  by  eddy  diffusion.  In  a  later  paper,  Anderson  . 
(1971b)  gives  a  measured  altitude  profile  for  hydroxyl  with  eryor  limits  1 
+  120  per  cent. 

At  present,  then,  estimates  of  chemical  abundances  of  minor  species 
required  to  predict  infrared  radiation  require  the  use  of  photochemical  • 
computations.  Most  computations  of  atmospheric  chemistry  before  the  past 
four  or  five  years  have  neglected  vertical  transport  by  molecular  diffusion 

and  eddy  mixing.  For  molecules  which  have  photochemical  lifetimes  greater 

A  5  i  !  : 

than  10  to  10  seconds,  this  neglect  may  yield  results  greatly  in  error,  as 

i 

a  comparison  of  the  calculations  of  Hunt  (1966),  Hesstvedt  (1968)  and  Shimazaki 
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mm 


and  Laird  (1970)  quickly  shows.  Even  calculations  including  transport 

■  .  :•  ' 

■  processes  cannot  be  accepted  at  face  value  because  the  values  of  eddy  mixing 

coefficients  used  at  present  are  for  the  most  part  only  informed  estimates, 

*  1.  • 

:  and  published  calculations  often  use  outdated  rate  coefficients.  Because  of 
the  above  reasons. and  the  future. need' to  estimate  abundances  for  conditions 
for  which  no  published  experimental  data  or  calculations  may  be  available, 
development  of  a  computational  model  for  atmospheric  photochemistry 

i 

including' vertical  transport  was  considered  desirable.  The  present  status 
of  this  model  is  described  below. 


i  Following  Shimazaki  and  Laird  (1970)  and  Bowman,  et  al.  (1970), 
the  continuity  equation  for  an  atmospheric  chemical  constituent  or  species 


may  be  written 


-istk  =  Q(z)  -  JUS)  -  ^ 
bt  az 


(2-1) 


toh^re  n  is  the  concentration  of  the  constituent,  Q  and  L  are  photochemical 
production  and  loss  rates  and  cf>  is  the  vertical  flux  due  to  transport  processes. 
When  only  vertical  transport  by  molecular  diffusion  and  eddy  or  turbulent 
mixing  is  considered,  the  vertical  flux  is 


’^>(z)  =  n(z)u(z)  +  n(z)v(z) 


(2-2) 


wherq  u  and  v  are  vertical  velocities  of  the  species  produced  by  molecular 
diffusion  and  eddy  mixing,  respectively.  If  D  is  an  '’average'*  molecular 
diffusion  coefficient,  H  is  the  scale  height  of  the  species  in  the  absence  of 


'  _ 


mixing  and  T  is  the  temperature, 


, .  _  an  ,  n  n  &T 

u(z  =  -  D  +  —  +  — ■  ~— 
iz  H  T  bz 


(2-3) 


If  K  is  the  eddy  mixing  coefficient  and  Hav  is  the  average  scale  height 
for  all  species  present. 


v<z)  =  -K  -is  +-a-+£41 

dz  H  T  d  z 
i  av 


(2-4) 


Then 


bn(z) 

dt 


^  =  01.)  -  L(.)  +  £  "f -gj  +  |§-  +5^-|  n  (2-i 


As  was  done  by  Bowman,  et  al.  (1970),  it  is  assumed  that  only  atomic 
oxygen  and  molecular  nitrogen  and  oxygen  are  important  as  collsion  partners. 
The  effective  diffusion  coefficients  for  atomic  oxygen  and  molecular  nitrogen 
and  oxygen  are  taken  to  be 


/k±- 

\  O  J  /  J  1J 


(2-6) 


where  Nq  is  number  density  at  standard  temperature  and  pressure  (STP), 
Tq  is  273.  15  K,  and  the  D„  are  mutual  diffusion  coefficients  at  STP.  The 


k^>sss;isi'®«35j55 


values  adopted  for  D^  „  ,  D_  _  and  D_  „  are  those  used  by 

°2‘N2  °"°2  °"N2 

Colegrove,  et  al.  (1966),  0.181,  0. 260  and  0. 260,  respectively.  That  for 
D^.  ^  ,  0.189,  is  taken  from  Chapman  and  Cowling  (1960),  while  Dx,  .. 

°2-°2  VN2 

and  Dq  q  are  estimated  to  be  0. 180  and  0.370.  The  effective  diffusion 
coefficients  for  other  species  are  taken  to  be 


D.  =  0.707  D_ 

i  o„ 


1  + 


m(Oz) 

m. 

i 


il/2 


(2-7) 


where  m.  represents  molecular  weight. 

In  converting  Equation  (2-5)  into  a  difference  scheme  for  numerical 
integration,  it  is  assumed  that  all  quantities  may  be  represented  by  a  quadratic 
function  of  altitude  at  any  given  altitude.  Then  ii  the  superscript  (+)  denotes  a 
value  at  the  next  higher  integration  point  and  the  superscript  (-)  denotes  a  value 
at  the  next  lower  point,  the  difference  equation  is 
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(2-8) 


The  computational  problem  then  is  simply  to  integrate  equation  (2-8) 
at  a  set  of  altitudes  with  spacing  Az»  taking  successive  time  steps  h  =  At 
over  the  desired  time  of  interest.  The  method  used  takes  advantage  of  the 
observation  that  except  where  a  species  is  important  as  a  third  body  in  an 
association  reaction  which  removes  it,  the  chemical  rate  equations  for  a 
coupled  system  may  be  cast  in  the  form 


x! 

1 


q.(x)  -^(xjx.  -  <*.(x)x. 


(2-9) 


with  formal  solutions  (assuming  q,  of  and  ^  constant  and  omitting 
obvious  subscripts) 


[(/3  +  S  )exp(-5t)  +  (<£-/0)]x.( 0)  +  2q[l  -  exp(-  &  t)3 
(2cCx.(0)  +/3)[l  -  exp(-  £  t)]  +  [1  +  exp(-  ^t)]^ 


q  >0 , y9  >0,  oC  >0 


(2-10a) 


5  = 


2  1/2 

(4c£q+/$  ) 


where 


where 


x{(t)  *  x{(0)txp(- ^t)  +  q[l  “  exp(-^t)]^# 


q>0f</5>0,  *.=  0  (2-lOh) 


[1  +  exp(-  £  t)]x.{0)£  +  2q[l  -  exp(-  £  t)] 
Xi(t)  =  2oCx.(0)[l  -  exp(-  & t)]  +  [1  +  exp(-£t)]5 


q>0,  /S=  0,  oC>0  (2-lOc) 


6=  (4oCq) 


x.(t)  -x.(0)  +  qt 


q>0,  y£=0,  <X=0 


/3x.(0)  expt-^t) 


(2-10d) 


Xi^  «<x.(0)[l  -  exp( 


-/*t)1  +  /3  ’  q  =  <K>  0  (2-lOe) 


x.(t)  =  x  (0)exp(-^$t),  q  =  0,  /£  >  0,  oC  =  0  (2-10*) 


x.(t)  =  x.(0)/[l  +  OCtx.(O)]  q  =  0,  /S  =  0,  «  >0  (2-10g) 


x.(t)  =  x.(0). 


q=  Of  /  =  0,  <*  =  0 


(2-lOh) 


10 


I 


Care  must  be  exercised  in  the  definition,  of  the  eC  and  ^  terms.  For 
example,  the  reaction  , 


O*  +  e  i  0  +  0  = 


i  i 

contributes  to  both  the  oC  and  the  /3  terms  in  the  differential  equation  for 

•f  ! 

Oz  concentration  if  other  positive  ions  are  present;-  > 


OC  =  k 


=  k([e]  -  [O*  ]) 


The  solutions.  Equation  (2-lp),  are- exact  only  V  and  oC  do  not 

change  during  an  integration  interval.  Since  q,  ^  ,  and  fC  are  functions  of  the 

time  dependent  concentrations  of  the  other  species  in  the  system,  the  actual 

"  .  *  ■ 

integration  scheme  must  allow  for  this  variation.  ’  i  1 

i 

If  the  time  step  t  2/ 1 2oc  x  +y9  -  £|  apd  ^t  «  2,  the  denominator 
of  Equation  (2- 10a)  may  be  expanded  as  1  -  (2 c<x  +/i-6)t/  26  and  , 


x.(t)  =  x.(0)  +  (q.  - /S.x.iO)  -  rf.x*  (0))  +  0(t2).  (2-11) 


For  a  small  enough  time  step,  the  algorithm  is  equivalent  to  Euler's  method. 

* 

With  the  same  restrictions  on  die  time  step,  the  algorithm  may  be  embedded  in 
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*■«*■-  ^>.v^,v>*  .-  £3$r<w3fl 


r~ 


-.  i 
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j 
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a  Runge-Kutta  formula  of  any  order  and  will  give  equivalent  results. 


.1 


■  The  program  uses  an  analog  of  the  improved  tangent  second-order 
Runge-Kutta  method  in  the  basic  step  h.  Let  x(o)  be  the  value  of  a  dependent 


M 

£ 


l£ 

>s 


variable  at  the  beginning  of  the  step.  Then  oC  (o),  A  (o)  and  q(o)  are  computed 

1  i  l '  . 

and  a  new  value  of  x  computed  for  t  =  -r-h  from  the  general  expression. 


Equation  (2-  10a),  or  from  an  equivalent  .expression,  if  any  oC,  ^3  or  q  is 

I  S 

equal  to  zero.  OC  (h/2),y£  (h/2)  and  q(h/2)  are  computed  and  Equation  (2-10) 

!  I  ,  «  » 

evaluated  with  these' values  and  t  =  h. 


Following  Ceschinotand  Kuhtzmann  (1966),  assume  that  the  stepsize  is 
chosen  so  that  a  linearized  theory  of  errors  may  be  applied.  That  is,  for  a 


method  of  order  n ,  the  error  for  one  step  of  length  h  is  Ah 


n+1 


Divide 


an  interval  Ht  into  an  even  number  k-  o(  steps  of  length  h .  Then 


.  (h).  (o)  ,.,n+l  ,  mun+2\ 

x.  (t)  =  x.  (t)  +  kAh  +-  0(h «  ) 


(2-12) 


1 


■  t 


where  x!  (t)  is  the  “exact”  result  of  the  integration  and  addition  or 


subtraction  of  the  error  by  propagation  is  neglected.  '  The  result  using  a  step 

<  S 

of  length  2h  is  1 


x!2h)(t)  =  x!0)(t)  +  k/2  A(2h)“'  +  Ofh*”  **)  (2-13) 


.n+1 


n+2. 


The  difference  between  the  :two  results  is 


12 


-5 

,3 


'll 

‘3 


$ 

I 


I 

I 


1 


-1 

I 


■sfcaiasM&fcgi 


x|2h)(t)  -  x|h)(t)  *  kAhn+1(2a-l) 


(2-14) 


and  the  principal  part  of  the  error  using  the  step  h  is 


£  *<2h)(0  -  *<»<t) 

«h=-i - 1 — 


2n  -  1 


(2-15) 


The  simplest  way  to  use  this  result  is  to  divide  an  interval  into  two 
equal  steps  h  and  integrate  twice,  once  with  two  steps  of  length  h  and  then 


with  one  step  of  length  2h .  The  principal  part  of  the  error  may  be  used  both 

_ 4.  _ -.14-  J2h) 


to  correct  the  result  x?  (t)  and  to  estimate  the  optimum  length  of  the  next 
interval  of  integration.  The  Runge-Kutta  integration  program  of  Gear  (1971) 
is  easy  to  adapt  for  both  uses  of  the  error  estimate,  and  was  used  in  early  tests. 


It  is  more  efficient  to  use  Romberg  extrapolation,  which  subdivides 
an  interval  H  in  more  than  two  ways,  and  uses  higher  order  error  terms. 
Either  polynomial  Romberg  extrapolation  [Davis  and  Rabinowitz,  (1967)]  or 
rational  Romberg  extrapol?.tion  [Bulirsch  and  Stoer  (1966)]  may  be  used. 


Denote  the  result  of  integrating  from  the  t  =  0  to  t  =  H  with  a  constant 
step  size  h  by  X?h).  Asstime  that  X(h)  has  an  asymptotic  expansion 


X(h)  =  5o  +  5jh 


V“ 


+  •  •  • 


(2-16) 


with  "S  ,  ...  constant,  so  that  in  the  limit  h  -*■(>,  X(0)  = 


13 


r 


the  exact  eolution.  An  approximation  to  X(0)  may  be  obtained  by  extrapolation 
from  the  reiulta  of  integration  with  finite  value#  of  h . 


Let  the  global  interval  of  integration  H  be  successively  divided  into 


2,  3,  4,  6,  8,  12,  ...  intervals  with  h  *  ^H,  ~H,  ^H,  ^H,  . . . 


Let  X*°^(H)  be  the  value  obtained  with  hs^H,  X^\h)  be  that 
o  J  Ct  O 

obtained  from  h  =  —  H,  etc.  Then  a  table  of  extrapolated  values 


X(o)(H) 

o 

X(1)(H) 

o 

X(j0)(H) 

X(2)(H) 

O 

X^tH) 

X^°;(H) 

X(3)(H) 

o 

X(2)(H) 

X^1}(H) 

X^0)(H) 

•  •  • 

X(m)(H) 

o 

«  •  • 

X1’n"1)(H) 

>4m‘2)(H) 

•  •  * 

(m-3) 

'3 

(H) 


is  computed  using 


r(o). 

kl 


X'  (H) 


A^X^H)  -  X^0)(H) 


a[o)  -  1 


14 


(2-17) 


..  X(t>)(H) 

1X1 


ffpilC &KKMB.M. _ 


>*» 


U>  A{1)X^2)(H)  -  X<1}(H) 

xi™  =  ^Tm - — 


,  .  A'  x'  (H)  -  X  (H) 
x'0)(H)  =  -2-  1  1 


-  1 


A 


(o) 


I 


(2-18) 


xj^(H) 


Ak>Xk-l)(H)  ‘  Xk-1(H) 


0) 


1 


If  the  method  used  in  integrating  over  a  step  h  is  of  order  n,  the 
exponents  in  Equation  (2-16)  have  the  values 


% .  =  n  -  1  +  i 
i 


(2-19) 


Equation  (2-16)  may  be  written  out  for  as  many  values  of  h  as  desired  and 

(j> 


these  equations  solved  simultaneously.  The  coefficients  Aj^  of  Equation  (2-18) 
are  obtained  from  the  multipliers  used  in  solving  such  a  system  of  algebraic 

(j> 


equations.  Table  II- 1  gives  the  coefficients  A£"  for  k  =  1,  2,  ...,  9  for 
first,  second,  third  and  fourth  order  integration  methods.  Only  the  second 
order  coefficients  are  used  here. 


Subdivision  of  the  basic  interval  is  continued  until  the  value  of 


15 


Table  II- 1 

Coefficients  A^  for  Polynomial  Extrapolation 
using  Subintervals  H/2,  H/3,  H/4,  H/6,  H/8,  ... 


First  Order 
i  even  j  odd 


Second  Order 
j  even  j  odd 


Third  Order 
i  even  j  odd 


Fourth  Order 


> 

O 

20 

14 

L* 

c 

7 

5 

% 

8 

13 

48 

J 

3 

3 

13 

A 

A 

40 

28 

T 

*X 

7 

5 

£ 

16 

198 

736 

D 

3 

23 

99 

8 

8 

80 

56 

7 

5 

12 

32 

292 

1088 

3 

17 

73 

16 

16 

160 

112 

7 

5 

24 

64 

3672 

13696 

3 

107 

459 

2320 

283 


3424 
20  9 

7716 

313 

43072 

1315 

3592 

73 


1132 

145 

3008 

291 


40064 

1929 

21040 

673 

18688 

449 


j  even 

j  odd 

8J. 

16 

256 

81 

208 

35 

70 

13 

173 

19 

i*»l4 

xoJ.u 

173 

1760 

149 

596 

55 

9666 

541 

69248 

4833 

23360 

991 

7928 

365 

89476 

2521 

645376 

22369 

10880 

231 

3696 

85 

11112 

80384 

16 


is  less  than  a  preassigned  value  €  .  In  practice,  adequate  results  are 
obtained  if  for  k  >6  the  test  is  on  values  of 


r(j> 


The  chief  advantage  of  this  integration  method  is  that  large  time 
intervals  may  be  used  and  the  algorithm  contains  an  error  estimate. 

As  applied  in  the  program,  the  L  term  in  Equation  (2-8)  is  divided 
2 

Into  C&y'.Hv  otfn  &r, „  To>  iuc  Q  c?rm  is  added  the  terms  containing 

•f  — 

n  and  n  .  The  multipliers  of  the  terms  in  n  are  added  to  the  y3  term 
from  the  photochemistry  and  Equation  (2-10)  is  used  with  interval  subdivision 
and  polynomial  Romberg  extrapolation  as  described  above. 

Table  II-2  lists  the  main  neutral  reactions  which  may  play  a  role 
in  the  chemistry  of  the  infrared  emitting  molecules  NO,  CO^,  H_,Ot  O^, 

N^O  and  CH^.  Table  II-3  lists  the  photodissociation  processes  considered. 
Table  II-4  lists  the  ionic  reactions  needed  for  nitric  oxide  chemistry.  Where 
possible,  an  estimated  error  is  given  for  the  value  of  the  rate  coefficient. 
Critical  reviews  are  the  principal  source  of  rates  listed  here  and  in  most 
cases  newer  rates  from  the  literature  fall  within  estimated  error  limits. 

Not  all  the  reactions  listed  in  Table  n-2  have  been  included  in  the 
program,  as  some  occur  in  branching  chains  which  have  not  been  verified  in 
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the  laboratory.  The  choice  of  reactions  ha:s  been  guided  by  those  tused  in 

i  I 

published  calculations  apd  reviews.  .  The  most  uncertain  rates  are  those 
involving  the  carbon  compounds.  These. have  been  studied  principally  at 
elevated  temperatures  and  extrapolation  to  the  low  temperatures  of  the  upper 
stratosphere  and  mesosphere  is  not  reliable.'  : 

The  results  of  Nicolet  (1970)  and  Bowman,  et  aTi^f  19.71)  were  used 
as  a  guide  to  the  selection  of  reactions  to  include  in  the  chemistry  of  oxygen 
and  hydrogen  compounds.  -  Rates  for  Reactions  Rl,  R2  and  R3  were  taken  - 

i 

from  Johnston  (1968)  in  preference  to  other  sources  because  they iare  thermo¬ 


dynamically  consistent  for  the  O-O^-O^  system.  The  rates  for  nitric  oxide 


reactions  follow  closely  those  used  by  Stroebel  (1971a,1 1971b).  The  chbice  of 
reactions  for  methane  decomposition  is  based  on  the  partial  kinetic  systems 

>  .  i 

of  Herron  (1969)  and  Hoare  and  Patel  (1969).  Reactions  leading  to  formation 


of  CH2  have  been  omitted. 


Of  the  chemical  species  considered  here,  nitric  oxide  abundances  are 


best  known  above  an  altitude  6f  70  km.  Observations  of  solar  ultraviolet 


induced  fluorescence  in  the  gamma  bands  have  been  made  by  Barth  (1964, 

»  •  i 

1966a,  1966b),  Pearce  (1969a)  and  Meira  (1971).  Abundances  measured  in 
this  way  are  much  higher  than  .predicted,  by  earlier  photochemical  theories 
such  as  those  of  Nicolet  (1965a;  1965b).  Norton  and  Barth  (1970)  and  Stroebel, 


2  ;  •  •  + 

et  al.  (1970)  showed  that  N(  D)  produced  in  dissociative  recombination  of  NO 


(Reaction  I14b)  and  in  charge  exchange  of  N_  with  atomic  oxygen  (Reaction  I7b) 

*•  i 

could  form  sufficient  nitric  oxide  in  reacting  with  molecular  oxygen  to  explain 


measured  abundances.  The  arhourit  of  N(  D)  formed  in  the  relevant  reactions 


has  not  been  determined  in  the  laboratbry  and  remains  an  adjustable  parameter 


in  all  computations  The  nitric  oxide  concentration  above  80  km  may  be 
considered  reasonably  well  known,  probably  to  within  a  factor!  of  two.  Below 


that  altitude,  the  chemistry  is  uncertain  and  the  effects  of  reactions  such  as 
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R54,  R55  and  R60  have  not  been  studied  thoroughly.  Stroebel  (1972)  shows 

*  1  i-  •  i  ,  ' 

that  values  lower  than  those  obtained  by  Meira  (1971)  arp  possible  and  suggests 

1  :  •  j  i  1 

using  limits  obtained  from  electron  densities  in  the  D  region.  Figure  II- 1 
shows  nitric  oxide  concentrations  used  here  as  a  basis  for  predicting  infrared 
radiance  at  high  altitudes,  and  is  based  on  Stroebel* s  calculations. 


Water  vapor  abundances  are  much  less  well  known.  ,There-is  general 
■  agreement  that  water  vapor  is  present  in, the  upper  stratosphere  in  amounts 
1  of  between  2  and  10  parts  per  million,  Anderson  (1971a)/  but  the  exact  amount 


is,  unknown.  Anderson  (1971a,  1971b)  has  determined  OH  abundances  from 


observation  of  fluorescent  scattering  of  solar  ultraviolet  radiation.  The, 


6  3 

concentrations  are  about  '5  x  10  molecules/cm  between  45  and  70  km,  with 


i  'i 

error  limits  of  about  120  per  cent.  Anderson's  first  paper  includes  a  steady  ; 
state  photochemical  model  with  diffusion  and  he  uses  this  to  infer  a  wafer  vapor 
’  abundance  ofabout  one  to  two  parts  per  million  at  85  km.  The  calculations 
with  the  chemistry  model  described  above  assumed  an  initial  concentration  of 
three  parts  per  million  and  were  not  run  long  enough  to  reach  diurnal  equilibrium. 
In  addition,  the  eddy  .diffusion  coefficients  above  50  km,  taken  from  Keneshea 

!  1  .  t 

and  Zimmerman  (1970),  are  larger  than  the  upper  limit  deduced  by  Anderson 
:(1971a).  It  was  decided  to  u se  them,  however,  because  M.  B.  McElroy 
(private  communication,  1972)  has  informed  us  that  there  has  been.aa 


unpublished  measurement  of  Reaction  R12  which  gives  a  value 


2  x  10  .^  cm^/sec  which  may  lead  to  lower  QH  concentrations  for  a  given 


water  vapor' abundaince.  The  adopted  water  vapor  abundances  are  given  in1 
Figure  fl-2,  and  are  based  oh  the  (incomplete)  chemistry  program  results. 

«  it  : 

i  * 

While  ozone  abundances  below  about  50  km'must  take  into  account 

*  *  ! 

lateral  transport  as  well  as  vertical,  transport  by  eddy  mixing  and  molecular 
diffusion,  photochemical  models  are  adequate  above  that  altitude,-  and  agree 


with  observations  obtained  by  measuring  absorption  of  solar  and  stellar 
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FIGURE  II- 1 


Altitude,  km 


ultraviolet  radiation  and  those  deduced  from  measurements  of  the  molecular 
oxygen  1.27  ^un  infrared  atmospheric  band.  The  exact  profile  depends  on 
assumed  eddy  diffusion  coefficients,  however,  principally  through  atomic 
hydrogen  and  oxygen  concentration.  abundances  used  here  are  shown  in 
Figure  II-3,  and  were  obtained  from  the  chemistry  program. 

Carbon  dioxide  has  a  lifetime  of  about  10  days  at  the  top  o:t  the  atmosphere 
where  it  is  dissociated  by  solar  radiation.  It  is  well  shielded  by  molecular 
oxygen  below  about  120  km,  however,  and  should  remain  well  mixed  below 
100  km.  The  only  good  measurement  above  100  km  is  that  of  Offerman 
and  von  Zahn  (1971),  who  find  that  between  120  and  140  km  the  mixing  ratio 
is  not  more  than  50  per  cent  less  than  the  tropospheric  value.  Above  140  km, 
molecular  diffusion  is  rapid  enough  compared  to  photodissociation  to  maintain 
diffusive  equilibrium.  While  Hays  and  Olivero  (1970)  have  obtained  a  rapid 
falloff  in  mixing  ratio  above  80  km,  their  results  are  incorrect  because  they 
neglect  recombination  by  Reactions  R33  and  R34  above  70  km.  While  the 
chemistry  program  was  not  run  long  enough  to  reach  realistic  values  of 
carbon  dioxide  concentrations,  the  profile  shown  in  Figure  II-4,  based  in  part 
on  Offerman  and  von  Zahn  (1971)  should  be  adequate  as  a  basis  for  computation 
of  infrared  radiances. 

Although  included  in  the  chemistry  program,  methane  is  not  considered 
as  an  important  source  of  infrared  radiation  above  50  km,  as  it  is  destroyed 
by  photodissociation  .  nd  reaction  with  atomic  oxygen  and  OH  faster  than  it  can 
be  transported  upward. 


Concentration,  molecules/cm 
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SECTION  m 


EXCITATION  MECHANISMS 


General  Considerations 


Except  for  the  pure  rotational  radiation  from  molecules  such  as 
water,  the  degree  of  excitation  of  vibrational  levels  determines  the 
radiation  from  infrared  emitting  molecules.  The  physical  processes 
affecting  the  population  of  vibrational  states  includes  collisional  excitation, 
radiative  excitation  by  resonant  absorption  or  fluorescent  decay,  chemical 
formation  in  a  vibrationally  excited  state,  and  radiative  and  collisional 


de-excitation. 


In  the  troposphere  and  lower  stratosphere  collisional  excitation  of 
vibrational  levels  is  sufficiently  rapid  to  compensate  for  radiative  de-excitation, 
and  vibrational  populations  are  in  thermodynamic  equilibrium  with  translation 
and  rotation.  At  altitudes  above  30  to  50  km,  however,  radiative  processes 
assume  a  dominant  role  in  determining  vibrational  populations.  Resonant 
scattering  and  fluorescent  mechanisms  become  the  dominant  processes 
producing  infrared  radiation.  At  wavelengths  greater  than  5 /cm,  the  most 


intense  source  of  radiation  is  the  thermal  emission  from  the  earth's  surface 


and  from  carbon  dioxide,  water  vapor  and  ozone  present  at  lower  levels  in 
the  atmosphere.  During  the  day,  the  effects  of  solar  radiation  must  be 
included  since  fluorescent  processes  may  lead  to  significant  excitation  of 
vibrational  states.  At  wavelengths  shorter  than  6  /cm,  resonant  absorption 
of  solar  infrared  flux  maintains  vibrational  temperatures  above  200  K. 


The  following  parts  of  this  section  discuss  vibrational  excitation  and 
de-excitation  and  die  model  developed  for  computation  of  radiative  processes. 


Preceding  page  blank 
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B.  Collisional  Excitation  ;  : 

In  computing  radiances  at  high  altitudes,  it  must  be  .ecognized  that 
vibrational  levels  are  not  populated  at  thermal  equilibrium  values.  Radiation 

de-populates  vibrational  states  rapidly  and  collisional  excitation  -is  not  . 

■ 

sufficiently  efficient  to  maintain  a  Boltzmann  distribution  of  vibrational 
states.  Collisional  excitation  and  de-excitation  rates  are  required  to  1 

■  j 

compute  the  balance  of  opposing  mechanisms  which  determine  actual 
vibrational  populations.  , 

Four  collisional  processes  are  considered  which  affect  the  vibrational 
state  of  a  molecule.  (1)  translational-vibrational  (T-V)  interactions  in  which 

1  * »  .  j 

tidnsiational  energy  is  transformed  into  vibrational  energy  or  vice-versa  during  , 
a  collision,  (2)  intermolecular  vibrational  exchange  (V-V)  in  which  a  molecule 

i 

of  one  species  exchanges  a  quantum  of  vibrational  .energy  with  a  molecule  of 
another  species,  (3)  intramolecular  vibrational  exchange  in  polyatomic 
molecules  in  which  a  collision  transfers  the  energy  in  one  mode  of  vibration  > 

c 

into  another  mode  within  the  same  molecule,  and  (4)  collisions  between  an  , 
atom  and  molecule  in  which  there  is  either  atom  exchange  or  the  formation  of 
a  collision  complex  which  has  a  lifetime  long  enough  fco  equilibrate  kinetic  and 
vibrational  temperatures.  This  latter  process  may  be  considered  a  special 
case  of  the  first  mechanism.  *  , 

I 

When  a  molecule  in  the1  ground  state  collides  with  another  molecule 
there  is  a  finite  probability  Pyj.  that  a- portion  of  the  relative  translational 
energy  will  be  transformed  to  vibrational  energy  if  the  translational  energy 
is  equal  to  or  exceeds  the  vibrational  energy.  Similarly,  if  the  molecule  is  in 
an  excited  vibrational  state,  there  is  a  finite  probability  P^,y  the  vibrational 
energy  will  be  converted  to  relative  translational  energy  during  the  collision. 

For  a  Boltzmann  distribution  of  translational  energies  at  a  given  temperature. 
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p  =  p 
TV  VT 
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(3-1) 


Sithilar  Expressions  hol’d  for  intermolecular  and  intramolecular  V-V  transfer. 

The  rate  at  which  a  vibrational  level  i3  excited  (k,^)  or  de_®xcited 
(kyT^  js  given  by  the  prqducfc  of  the  collision  frequency  Z  and  the  transition 
probability.  Experimental  determinations  of  collisional  excitation  rates  may 
be  reported  as  probabilities,  excitation  r^tes  or  relaxation  times.  The  latter 
are  related  to  excitation  and  de-excitation  rates  by  (Herzfeld  and  Litovitz,  1959) 


—  =  k„_  +  k 
T  i  VT  TV 


(3-2) 


Excitation  rates  depend  on  temperature  and  number  density  through  the  strong 
temperature  variation  of  the  transition  probability  and  the  square  root  tempera¬ 
ture  and  linear  density  variation  of  the  collision  frequency. 

•  i 

The  most  important  molecular  collision  partners  in  the  atmosphere  are 
« 

nitrogen  and  oxygen,  and  a  knowledge  of  their  vibrational  temperatures  and 
effects  is  important.  These  will  be  treated  first  and  followed  by  a  brief 
discussion  of  the  relative  importance  of  collisional  excitation  mechanisms 
for  nitric  oxide. 

i 

Vibrationally  excited  oxygen  can  be  produced  in  the  atmosphere  by 
recombination  of  atomic  oxygen 


0  +  0+  M— ►Oj  +  M 


(3-3) 
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by  formation  from  ozone 


0  +  0, 


°2  +  °2 


by  fluorescence  in  the  discrete  Schumann-Runge  band  system 


hV  +  O^  -*• 


°2  +  h* 


by  transfer  of  kinetic  energy  to  vibrational  energy 


Oz  +  M  -►  Oz  +  M 


by  electron  impact 


e  +  oz  —  e  +  02 


by  exchange  of  vibrational  quanta  with  molecular  nitrogen  (and  other 
molecules) 


N2  +  °2  -  NZ  +  °2 


and  by  the  strong  collisionai  interaction  with  atomic  oxygen 


O  +  o_  0_  +  O  . 


VibrationaUy  excited  oxygen  can  be  quenched  by  the  strong  collisional 
interaction  with  atomic  oxygen 

O  +  O*  — -  O  +  Oz  (3-10) 

by  transfer  of  vibrational  quanta  to  molecular  nitrogen  (and  other  molecules) 

°2  +  N2  —  °2  +  N2  <3-U> 
and  by  transfer  of  vibrational  energy  to  kinetic  energy 

O*  +  M  Oz  4  M  (3-12) 

Above  the  mesopause,  and  up  to  300  km,  reactions  (3-9)  and  (3-10) 
dominate  production  and  loss  of  vibrational  excitation  and  can  be  shown  to 
maintain  vibrational  distributions  close  to  thermal  equilibrium  values. 

Bates  and  Moiseiwitsch  (1956)  proposed  that  the  atom  exchange 
reaction 


O'  4-  Oz  —  O'O  +  O  (3-13) 

is  sufficiently  rapid  in  quenching  vibrationally  excited  oxygen  that  the  reaction 
proposed  by  Krassovsky  (1955) 

O*  4-  H  — *■  OH* 4-  O  (3-14) 

cannot  be  important  in  the  production  of  vibrationally  excited  OH.  Bauer  and 
Tsang  (1963)  suggested  that  reaction  (3-14)  would  be  important  in  relaxing  the 
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vibrational  levels  of  shock  heated  oxygen  and  estimated  the  rate  coefficient 

-10  2 
to  be  about  2  x  1C  exp  (- 8000/RT)  cm  /sec. 

The  rate  for  atom  exchange  reaction  (3-14)  has  been  measured  near 

room  temperature  by  Ogg  and  Sutphen  (1953,  1954),  Brennan  and  Niki  (1965), 

Jaffe  and  Klein  (1966)  and  Herron  and  Klein  (1964).  Johnston  (1968)  has 

-12  3 

reviewed  these  results  and  recommends  a  value  of  10  cm  /sec  for  the 
rate  coefficient  at  273  K.  Garnett,  Kistiakowski  and  O'Grady  (1969)  have 
measured  the  rate  coefficient  near  1800  K.  Clark,  Garnett  and  Kistiakowski 
(1970)  have  combined  the  latter  results  with  those  of  Jaffe  and  Klein  (1966)  to 
suggest  a  temperature  dependence  given  by 

i:i3  =  10'12*48-°,3T1/2exp[-(960  +  700)/RT3Cm  /s6C  (3-15) 

Kiefer  and  L-utz  (1967)  have  measured  the  vibrational  relaxation  time 

of  molecular  oxygen  in  the  presence  of  atomic  oxygen  at  temperatures  between 

1600  and  3300  K.  They  find  that  over  this  temperature  range  the  product  of 

pressure  P  and  relaxation  time  "£  is  given  by 

°2 

FT  =  (4.35  +  0.  19)  x  10"8  -  (7.75+0.81)  x  10"12T  atm-sec  (3-16) 
2 

The  prediction  of  Bauer  and  Tsang  (1963)  agress  with  this  result  within  a  factor 

of  two.  The  rate  coefficient  for  atom  exchange  from  Equation  (3-15)  is 
- 12  3 

9.8  x  10  cm  /sec  at  1600  K  and  the  de-excitation  rate  coefficient  derived 

-12  3 

from  Equation  (3-16)  is  8.6  x  10  cm  /sec  at  1600  K.  The  close  agreement 


mmmmmxemmswsma 
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g  =  - 7  TTF  f  A 


(3-17) 


supports  the  assumption  used  here  that  the  vibrational  de-excitation  rate 
is  within  a  factor  of  two  of  the  atom  exchange  rate,  which  may  be  considered 
well  determined. 

Further  evidence  is  given  by  the  theoretical  computation  of  Breig 
(1969)  who  estimated  the  rate  of  de-excitation  of  high  vibrational  levels  of 
molecular  oxygen  by  reaction  with  atomic  oxygen.  He  obtains  total  reaction 
rates  of  the  V  =  5  and  V  =  25  levels  at  200  K.  Breig*  s  (1969)  results 
may  be  used  as  a  qualitative  argument  for  the  further  assumption  that  higher 
vibrational  levels  of  molecular  oxygen  will  reach  steady  state  populations  in 
times  of  the  order  of  that  required  by  the  V  =  1  level. 

Oalgarno  (1963,  1967,  1970)  has  suggested  that  fluorescence  in  the 

discrete  Schumann-Runge  band  system  may  be  an  important  source  of 

vibrationally  excited  oxygen  in  the  mesosphere  and  thermosphere.  The 

experimental  results  of  Ackerman,  Biaume  and  Nicolet  (1969)  and  Hudson  and 

Carter  (1969)  indicate  that  excitation  to  V*>  3  results  in  predissociation 

rather  than  fluorescence.  If  this  conclusion  is  accepted,  the  only  transition 

o 

of  importance  is  the  (2,  0)  transition  at  1971  A.  This  has  an  oscillator  strength 

-4 

f  in  absorption  of  2.34  x  10  according  to  the  computations  of  Nicnolls  (1964). 

o  I]  2  o 

The  solar  flux  "TTF  at  1971  A  is  about  1.2  x  10  photons/(cm  -sec-A) 

°  11  2  ° 
according  to  Detwiler,  et  al.  (1961)  or  1.4  x  10  photons/(cm  -sec-A) 

according  to  Ackerman  (1970).  At  the  top  of  the  atmosphere,  the  absorption 

_7 

rage  g  per  oxygen  molecule  in  the  (2,  0)  band  is  9.6  x  10  photons  per  second 
from  the  relation 


Using  the  Einstein  A  coefficients  given  by  Nicholls  (1964),  it  is-found  that 

i 

the  mean  vibrational  level  following  fluorescent  emission  is  11  so  that  about 

-5  -  . 

10  vibrational  quanta  are  produced  each  second  per  molecule. 

i 

The  effective  rates  of  vibrational  excitation  by  electron  impact  may 
be  determined  from  the  electron  cooling  rates  computed  by  Lane  and 
Dalgarno  (1969).  Using  the  electron  densities  and  temperatures  observed 
by  Spencer,  et  al.  (1965),  ^typical  rates'are  3.  3  x  10  T,  8  x  10  and 
6.2  x  10  ^  excitations  per  molecule  per  second  at  150,  200  and  3001  km, 
respectively.  : 

The  effects  of  the  above  processes  on  molecular,  oxygen  vibrational 
temperatures  were  computed  for  the  CIRA  1965  mean  atmosphere  for  altitudes 
between  80  and  300  km.  The  rate  coefficients  given  by  Johnston  (1968)  were 
used  for  three  body  recombination  of  atomic  oxygen  and  for  the  reaction  of 
atomic  oxygen  with  ozone.  The  jnaximum  possible  vibrational  excitation  of 
the  products  was  assumed.  The  rate  coefficients  given  by  Gilmore,  Bauer  ' 

and  McGowan  (1969)  were  used  for  vibrational  excitation  by  collisional: 

! 

transformation  of  kinetic  energy  and  collisional  exchange  of  vibrational  quanta 

I  1 

with  molecular  nitrogen.  The  de-excitation  rate  by  reaction  with  atoiqic 

•  _i3  1/2  3 

oxygen  was  taken  to  be  3.  3  x  10  |  T  exp(-485/T)  cm  /sec,  the  mean 
value  suggested  by  Clark,  Garnett  and  Kistiakowski  (1979). 

:  i 

Even  with  the  upper  limits  to  nitrogen  vibrational  temperatures  given 
by  Walker,.  Stolarski  and  T'lagy  (1969)  thft  rate  of  vibrqtional  excitation  of 
molecular  oxygen  by  collisional  transfer  of  energy  with ‘nitrogen  is  less  than 
19  ^  excitations  per  molecule  per  second  ‘at  any  altitude  above  80  km. 

i 

1 

At  all  altitudes  below  200  km,'  die  vibrational  population  of  molecular 

*  1 

oxygen  as  computed  is  within  two  per  cent  of  the  thermal  equilibrium  value. 

i  . 

At  300  km,  excitation  by  electron  impact  increases  the  population  of  the  first 
vibrational  level  by  less  than  twenty  per  cfeht.  Thus,  at  all  altitudes  where 


.  1  '  !  ‘  '  "* 


polar  molecules  will  be  significantly  excited  by  molecular  oxygen,  the  ! 
vibrational  temperature  may  be  taken  to  be  the  kinetic  temperature.  This 

i  s 

resblt  agrees  with  the  conclusions  df,  Bauer,  Kummler  and  Bortner  (1971). 

*  i 

The  vibrational  temperature  of  molecular  nitrogen  at  high  altitudes 

I  :  j  } 

remains  uncertain.  The  principal  processes  which  have  been  suggested. to 
■  •  1  :  ‘  ' 

produce  vibrationally  excited  nitrogen  are  transfer  of  kinetic  energy  to  . 

•  i 

vibrational  auanta 


N2  +  M 


N2  +  M 


(3-18V 


quenching  of  0(*D)  by  collision 


Of1©}  +  N2  -♦  0(3P)  +  N2(V£7)  '  (3-19} 


reaction  of  atomic  nitrogen  and  nitric  oxide 


N(43)  +  NO  N*  (V  <  12}  +  0(3P) 

‘  w 


(3-20) 


and  excitation  by  electron  impact 


e  +  N. 


$ 

e  t  *2 


(3-21  )' 


The  principal  propess  which  quenches  vibrationally  excited  nitrogen  is 
collisional  transfer  of  vibrational  quanta  to  the  vibrational  mode  of 
carbon  dioxide 


N*  +  COz  —T  N,  +  :C02 


(3-22) 
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Reaction  (3-18)  is  insignificant  above  70  or  80  km.  If  all  the 
available  energy  were  transformed  to  vibrational  excitation,  quenching  of 
0(*D)  would  be  the  most  important  process  above  80  km.  However,  Fisher  and 
Bauer  (1970)  estimate  that  the  process  is  only  about  five  per  cent  efficient  and 
Black  and  Eckstrom  (1971)  have  obtained  an  experimental  efficiency  of 
8.  3^6.5  per  cent.  Above  120  km,  reaction  (3-20)  is  the  most  important 
source  of  excitation,  Morgan,  et  al.  (1962)  and  Morgan  and  Schiff  (1963) 
having  determined  that  about  one-third  cf  the  available  energy  goes  into 
vibrational  excitation.  Reaction  (3-21)  is  nearly  as  an  important  source, 
the  electrons  involved  being  those  produced  by  photo-ionization  in  the 
E- region. 

Walker  (1968)  and  Walker,  Stolarski  and  Nagy  (1969)  have  computed 
theoretical  estimates  of  nitrogen  vibrational  temperatures.  The  first  paper 
neglected  collisional  transfer  of  vibrational  energy  to  carbon  dioxide  and 
obtained  an  approximately  constant  vibrational  temperature  of  3000  K  above 
110  km.  Corbin,  et  al.  (1970)  noted  that  this  result  was  too  high  near  110  km, 
more  energy  being  transferred  to  carbon  dioxide  than  could  be  produced  by 
the  assumed  excitation  mechanisms.  The  second  paper  included  quenching 
by  carbon  dioxide  and  recognized  that  the  0(*D)  excitation  mechanism  might 
be  inefficient.  Results  were  given  for  zero  and  100  per  cent  excitation  by 
this  mechanism.  The  first  assumption  gives  a  maximum  vibrational  temperature 
of  1700  K  at  300  km.  The  vibrational  temperature  is  below  1000  K  at  all  times 
for  altitudes  below  140  km.  Minimum  vibrational  temperatures  occur  at 
dawn  and  rise  steadily  until  sunset.  Using  these  results,  it  is  found  that 
excapt  for  carbon  dioxide  band  emissions  originating  from  the  mode, 
excess  vibrational  excitation  of  nitrogen  is  unimportant  for  radiation  emitted 
by  the  molecules  studied  here. 

Degges  (1971)  has  discussed  the  mechanisms  which  contribute  to 


vibrational  excitation  of  nitric  oxide.  He  proposes  that  above  80  km, 
the  most  important  collisional  process  is  the  atom  exchange  mechanism 


NO  +  O' 


NO*  +  O 


(3-23) 


Using  the  low  temperature  measurement  of  Herron  and  Klein  (1964)  and  an 

argument  based  on  the  adiabatic  correlation  diagram  of  Donovan  and 

Husain  (1970)  and  the  temperature  dependence  of  the  similar  atom  exchange 

reaction  of  atomic  oxygen  with  carbon  monoxide,  he  suggests  that  the 

de- excitation  rate  of  vibrationally  excited  nitric  oxide  by  collision  with 

-13  1/2  3  -13  1/2 

atomic  oxygen  lies  between  10  T  cm  /sec  and  5  x  10  T  exp(-1000/RT) 

3 

cm  /sec.  The  latter  race  has  been  adopted  in  the  present  computations. 

Figure  HI-  lcompares  the  effect  of  this  process  with  exchange  of  vibrational 
quanta  with  molecular  oxygen  and  nitrogen.  Even  if  the  nitrogen  vibrational 
temperature  is  2000  K,  nitrogen  is  a  less  efficient  collision  partner  above 
120  km. 

The  rates  used  for  other  vibrational  excitation  processes  are  largely 
the  same  as  those  used  by  Corbin,  et  al.  (1969).  Vibrational  exchange  rates 
between  nitric  oxide  and  nitrogen  are  those  measured  by  Basco,  et  al.  (1961, 1962) 
with  the  temperature  dependence  given  by  Gilmore,  et  al.  (1969)  which  is  based 
on  the  correlation  scheme  of  Millikan  and  White  (1963).  Collisional  excitation 
rates  used  for  carbon  dioxide  and  water  vapor  are  taken  from  the  review  of 
Taylor  and  Bitterman  (1969),  with  the  addition  of  the  experimental  results  of 
Stephenson  et  al.  (1971)  on  intramolecular  vibrational  exchange  in  carbon  dioxide. 
The  rates  used  for  nitrous  oxide  are  those  of  Yardley  (1968)  and  Simpson,  et  al. 
(1968).  Estimated  rates  for  ozone  excitation  were  obtained  by  use  of  Millikan  and 
White* 8  (1963)  correlation  scheme,  and  comparison  with  other  systems.  Because 
of  the  lower  temperatures  at  the  mesopause  (about  80  to  90  km)  and  the  rapid 
decrease  in  number  density  with  altitude,  the  only  critical  rates  above  80  km  are 
the  oxygen  atom  exchange  with  nitric  oxide  and  the  vibrational  exchange  between 
carbon  dioxide  and  molecular  nitrogen. 
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Radid  i  ive  Processes  and  Modelin; 


The  physical  processes  and  assumptions  used  in  developing  a  model 
for  transport  of  infrared  radiation  in  the  upper  atmosphere  are  discussed 
below,  beginning  with  the  interaction  of  electromagnetic  radiation  with  an 
isolated  atom  or  molecule.  An  understanding  of  the  basic  processes  involved 
is  oie  first  step  in  any  computation  of  the  radiation  from  a  gas  or  mixture  of 
gases.  Subsequently,  absorption  from  a  continuum  source,  a  spectrally 
varying  source  and  fluorescent  processes  is  treated.  The  effects  of  optical 
thickness  when  the  path  of  the  radiation  includes  a  large  number  of  molecules 
of  the  same  species  is  considered  for  spectral  lines  with  Doppler  shape. 
Finally,  a  model  for  the  transport  of  radiation  for  a  single  vibration-rotation 
band  is  described. 

When  treating  a  problem  in  which  the  number  of  photon  emitted  or 
absorbed  by  an  atom  or  molecule  must  be  considered  rather  than  total  energy 
emitted,  it  is  convenient  to  examine  the  physics  of  the  interaction  of  radiation 
with  matter  through  use  of  the  Einstein  A  and  B  coefficients.  In  treating  a 
radiative  problem  such  as  the  current  study  of  high  altitude  infrared  radiances 
which  involves  transfer  of  energy  between  molecules  through  collisions,  this 
approach  is  a  necessity.  In  the  following  discussion,  cgs  units  are  used 
except  where  explicitly  stated. 

The  Einstein  A^  coefficient  is  the  rate  at  which  an  atom  or  molecule 
in  an  excited  state  spontaneously  emits  radiation.  It  has  the  units 
photons/ sec-molecule.  In  an  assemblage  of  molecules  all  in  a  state  given  by 
the  Einstein  A^j  coefficient  for  the  particular  radiative  transition  from  a  state 
of  higher  energy  (subscript  u)  to  one  of  lower  energy  (subscript  1). 

In  addition  to  spontaneous  emission,  two  other  processes  occur-  These 
are  absorption,  characterized  by  the  Einstein  coefficient  for  '.bsorption  5,  , 


*rt£VU4«t  Wfeh«UMM0M-rtM«MWrW  »*••••*< 


and  induced  emission,  characterized  by  the  Einstein  coefficient  for  induced 


emission  B  .  .  In  a  radiation  field  with  spectral  density,^^  ,  the 
•  111  lu 


probability  that  an  atom  or  molecule  absorbs  a  quantum  of  radiation  in  unit 


time  is  B,  ^  •  When  radiation  density  is  expressed  in  c.g.  s.  units 

lu’  -V  3 

ergs-sec/cm^,  coefficient  B^  has  units  cm  /ergs-sec.  Similarly, 


an  atom  or  molecule  in  a  radiation  field  with  the  above  spectral  density  and 

in  an  excited  state  emits  radiation  at  the  rate  E^/^y  *n  Edition  to  that 

lu 

characterized  by  the  coefficient  A^.  The  Einstein  coefficients  are  properties 
of  an  atom  or  molecule.  If  one  coefficient  is  known,  the  others  may  be 
determined  by  the  relations 


A  .  =  81ThB  JX" 
ul  ul 


(3-24) 


and 


g.B  =  s.  B  i. 

1  lu  u  u 


(3-25) 


Thus,  it  is  possible  to  treat  the  problem  knowing  only  the  Einstein 
coefficient. 


Experimental  determinations  of  the  interaction  of  infrared  radiation 


-:5t 


i 


with  gases  are  usually  made  by  measurements  of  light  absorption.  Experimental  >  ® 

'  H 

results  may  be  expressed  as  in  ter?  I'd  ted  absorption  from  a  continuum  source,  ’’f¬ 


lu 


=  N. 


A  l2  r 

A  .A.  g 
ul  lu  u 


1  81T 


g 


1 


1  -  exp 


(3-26) 
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or  as  line  or  band  strengths 


S,o  •  Sia/pc  (3-27) 

In  the  visible  and  ultraviolet  regions  of  the  spectrum,  it  is  customary 
to  use  oscillator  strengths  to  characterize  the  interaction  of  radiation  and 
matter.  In  this  case,  the  defining  equation  is 


s;  =  ~ -  N  f 

lu  me  1  Ju 


1  -  exp 

he  \1 

['  A,  kT 

l  lu  /  J 

(3-28) 


The  results  of  theoretical  computations  of  the  interaction  of  matter 
with  radiation  are  frequently  expressed  as  dipole  moments,  These 

are  related  to  the  Einstein  A  coefficient  by 


=  64  it4  h 

3hA3  gu 


(3-29) 


The  rate  at  which  a  single  atom  or  molecule  absorbs  and  reradiates 
light  from  a  parallel  beam  from  a  continuum  source  at  a  wavelength  Aq 
corresponding  to  a  resonance  transition  is  (Mitchell  and  Zemansky,  l$'6l): 


ITe2  2 

2  Vo  *o 

(photons/sec) 

me 

8.852  x  10_13I  f  A2 

n  n  n 

(photons/sec) 

(3-30) 
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where  I  is  the  continuum  source  flux  alnd  f  is  the;  oscillator  strength. 


If  the  incident  flux  is  expressed  in  photons/cm2- jjlxti,  and  oscillator  strength 

•  -2-1 
is  replaced  by  the  equivalent  line  or  band  strength  S  expressed  as  cm  ..  atm  ,, 


1 


G  =  3.72  x  10 "24I  SA2 

o  : 


i  (3-31)’ 
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Implicit  in  this  expression  is  the  convenient  relation  between  the  Einstein 
coefficient  and  line  or  band  strength  with  wavelength  ’dimensions  juurx 


=  2.804  S{cm"  2atm-1)/ A  2(>U-m).  (3-32) 


-1 


When  the  source  is  an  extended  black  body  surface  with  a  spectral 


radiance  ,  the  G  factor  is  obt*  ined  by  intergration  over  the  solid  angle 
subtended  by  the  source.  For  radiation  from  a  plane  parallel  surface  of 
infinite  extent. 


G  =  1.169  *  10"  SA 


(3-33) 


Another  purely  radiative  mechanism  of  importance  in  upper  atmosphere 

i 

radiance  computations  is  fluorescence.  A  striking  example  occurs  for  the 
weak  10.4  fj. m  CO.>  band.  Absorption  of  solar  radiation  by  CO.,  4.3  yttm  band 


V3  vibrational  mode  at  a-  vibrational  temperature  of  about 


maintains  the 

^  ! 

290  K.  During  the  day,  this  leads  to  an  enhancement  of  emission  in  the 


'4* 


■Passim 


imasa  sis  .. 


10.4  fun  band  by  as  much  as  two  orders  of  magnitude  when  viewed  along  the 
earth's  limb- from  outside  the  atmosphere. 

Fluorescence  ocdurs  when  the  upper  state  of  the  resonance  transition 
reradiates  into  other  than  the  initial  states.  The  relative  probability  of  a 
radiative. transition  to  a  lower  state,  n  from  the  upper  state,  u  when  there 

i 

are  k  possible  lower  states  is: 

i 


i 


(3-34) 


where  the  A  .  are  the  Einstein  A  coefficients  for  the  spontaneous  emission 
uj  i 

probabilities  for  a  transition  from  the  upper  state  to  the  jth  lower  state. 

1 

In  the1  case"  of  molecular  resonance-fluorescent  scattering  at  low 

.  » 

temperature,  only  the  lowest  vibrational  level  (v"  =  0)  in  the  ground  electronic 
state  is.  populated  significantly.  Absorption  of  a  photon  excites  the  molecule 
to  a  vibrational  lev4l,  v*  in  the  same,  or  a  higher  electronic  state.  Radiation 
is  possible  from  this  state  to  the  initial  level  (v"  =  0)  or  to  other  vibrational 
levels  (v"  ii  0)  in  the  ground  electronic  state.  The  rate  of  resonance  fluorescent 
scattering  in  a  given  electronic  vibrational  transition  may  be  obtained  by 
combining  Equations  (3-30)  and  (3-34)  and  introducing  appropriate  subscripts: 


v'v'f 


8.853  x  10-1?I 


v*  v" 


v'o  v'o 


Av'v" 


(3-35) 
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Here  the  subscript  v'v11  refers  -o  a  transition  from  an  upper  state  with 
vibrational  level  v’  to  a  lower  state  with  vibrational  level  v'*.  The  subscript 
v*o  refers  to  absorption  of  a  photon  by  a  molecule  in  the  lowest  vibrational 
le^eJ.  of  the  ground  state,  exciting  the  molecule  to  a  vibrational  level  v*  in 
the  same  or  a  higher  electronic  state. 

In  developing  a  realistic  model  of  the  interaction  of  radiation  with 
matter,  it  is  necessary  to  tak-*  into  account  the  details  of  absorption.  Line 
shape  must  be  considered  because  it  determines  the  amount  of  radiation 
escaping  from  an  optically  thick  layer  of  gas.  Thus,  for  large  optical 
thicknesses  (greater  than  a  magnitude  of  about  10),  the  total  emission  from 
an  isolated  spectral  line  with  Lorentz  shape  is  proportional  to  the  square 
root  of  the  optical  thickness,  while  that  from  a  line  with  Doppler  shape  is 
proportional  to  the  square  root  cf  the  logarithm  of  the  optical  thickness.  In 
addition,  absorption  of  radiation  emitted  at  one  place  in  the  gas  by  a  molecule 
at  another  location  cannot  be  computed  without  knowledge  of  the  line  shape. 

At  high  altitudes,  the  profile  of  the  spectral  absorption  curve  of  an 
individual  line  in  a  vibration-rotation  band  is  predominantly  due  to  velocity 
or  Doppler  broadening.  Kuhn  and  London  (1969),  who  investigated  the  heat 
budget  of  the  atmosphere  between  30  and  110  km,  report  that  radiative 
transfer  calculations  require  the  use  of  the  Voigt  profile  (including  both 
collisional  and  Doppler  broadening)  at  altitudes  between  20  and  70  km. 
Alternatively,  the  Doppler  profile  is  adequate  at  higher  altitudes  where  the 
Doppler  shape  is  the  largest  contributor  to  line  broadening.  The  basic 
computational  model  therefore  assumes  a  Doppler  profile. 

In  addition  to  the  use  of  Doppler  shape  for  a  single  line,  several 
other  assumptions  have  been  made  about  the  radiative  transfer  within  lines 
of  a  band.  These  are  the  following: 
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(1) 


Rotational  levels  within  a  vibrational  level  remain  in 
equilibrium  with  the  translational  temperature.  This 
is  valid  up  to  at  least  150  km  according  to  Goody  (1964). 


(2) 


There  is  complete  redistribution  of  frequency  within  a 
single  line  following  absorption  of  radiation. 


(3) 


There  is  complete  redistribution  of  absorbed  energy 
among  all  rotational  levels. 


(4) 


Line  strengths  are  used  as  if  each  line  were  at  the 
band  center,  neglecting  wavelength  variation  across 
a  band. 


(5) 

(6) 


There  is  no  overlap  of  adjacent  lines. 


Temperature  variations  at  different  levels  do  not  affect 
the  radiative  transfer.  This  assumption  probably 
constitutes  the  largest  source  of  error  in  computing 
transmission. 


The  major  result  is  to  make  computed  band  radiances  too  small  when 
determining  escape  of  radiation  from  a  warm  layer  overlaid  by  a  cooler 
layer. 


The  above  assumptions  lead  to  the  following  treatment,  which  begins 
with  consideration  of  a  single  line  and  from  this  develops  the  method  used  for 
treating  radiative  transfer  in  an  entire  band.  The  spectral  absorption 
coefficient  of  a  Doppler  broadened  line  is  (Mitchell  and  Zemansky,  1961): 


k(V)  = 


Avt 


In  2  tfe,T  t 

- — - N„f  exp 

me  1  o 
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fin 


(3-36) 
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and  the  total  light  absorbed  is 


r«>  r°° 

j  I(^)A(V)dV  =  I(y)  (1  -  T(>>))  d>> 
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The  quantity  in  parentheses  in  Equation  (3-40)  is  simply  a  factor  .by 
which  the  product  of  G  and  must  be  multiplied  to  account  for  the  effects’ 
of  total  optical  thickness  in  computing  the  total  absorption  of  radiation 

by  a  column  containing  Nj  molecules/cm2  in  the  lower  state.  This  quantity 
is  termed  S  by  Mitchell  and  Zemapsky  (1961)  who  present  tabulations  of  its 

■  i  ■ 

value.  Penner  (1959)  gives  an  asymptotic  expansion,  for  large  values  of  fC^. 
The  astrophysical  literature  frequently  uses  the  notation 


L(  K  ) 

o 


■v  .-4T- 


(3r41) 


The  factor  by  which  G  must  be  multiplied  to  obtain  the  true  rate  of 

i  *  i  • 

absorption  at  an  optical  depth  may  be  found  by  differentiating  Equation 
(3-41)  with  respect  to  K  (Ivanov  and  Shcherbakov,  1965): 

O  !* 


dL(  KJ 

Mi<  V  =  ~nr 

o 


’2  -w2  -  K^e-W 

w  e 


(3-4Z) 


The  factor  by'which  G  must  be  multiplied  to  obtain  the  actual 
absorption  from  an  extended  black  body  source  is  (Ivanov  and  Shcherbakov, 

t 

1965) 


2tT  JV2 


r  f 

W  =  (  ** 


cos  6  sin©  M 


11  cos  B , 


(3-43) 
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When  the  source  of  radiation  is  a  Doppler  profile  exactly  that  of 
the  absorbing  line,  the  differential  absorption  at  optical  thickness  H 


is  given  by 


s*  V  ir[  e  ’ 


(3-44) 


and  G  must  be  multiplied  by  this  fabtor.  Also,  if  t!be  source  is  an  infinite 
plane  layer  (not  a  surface)  radiating  with  a  Doppler  profile,  G  must  be 

i  ; 

multiplied  by* 


N2(  *o>  *  [ 


2  If  M  2 


,Vcos«  / 


sin0  dO 


(S'45) 


Corbin,  et  al.  (1970)  evaluated  the  functions  of  Equations  (3-41)  , 

I  ' 

through  (3-45)  by  numerical  intergration  and  obtained  polynomial  ■ 
approximations  for  the  entire,  range  of  variation  of  xq  from  zero  to 
infinity,  usiug  the  known  asymptotic  forms  as  a  guide  at  large  values.  They 
then  used  these  functions  to  compute  the  corresponding  functions  for  vibration- 
rotation  bands  of  linear  molecules.  The  resulting  functions  are  conveniently 

*  s  » 

parameterized  by  using  the  optical  thickness  at  the  center  of  the  strongest 
line  of  a  band.  It  was  found  that  two  sets  of  such  functions  were  adequate, 

!  !  i 

one  set  for  parallel  transitions  and  one  set  for  perpendicular  transitions. 

.  I 

» 

In  extending  the  modeling  to  non-linear  molecules,  it  was  found  to  be 

.  *  '  *i  * 

more  convenient  to  use  the  number  of  molecules  in  the  line  of  sight  as  the 


:(^t^533i£=5s-5Wf‘.=h 


model  parameter.  Since  upper  atmosphere  temperatures  are  near  or 
below  250  K  where  optical  thickness  affects  results,  the  functions  used  in 
the  program  for  the  water  vapor  ozone  Vj  and  bands  were 

computer  for  a  temperature  of  250  K.  Because  ozone  and  water  vapor  are 
present  in  much  smaller  concentrations  than  carbon  dioxide,  the  effects  of 
optical  thickness  are  mvch  smaller. 


At  altitudes  where  collisional  excitation  rates  are  small,  particularly 
between  80  and  100  km  and  above  120  km,  the  radiation  from  the  lower 
atmosphere  and  the  sun  will  control  vibrational  excitation  of  infrared  emitting 
molecules.  Figure  III-2  shows  vibrational  temperatures  which  will  be  maintained 
by  radiation  from  the  earth's  surface  and  lower  atmosphere  assuming  that  the 
upward  radiation  may  be  treated  as  originating  from  a  blackbody  at  a  given 
temperature.  The  effective  cemperature  will  vary  with  wavelength  and  lower 
atmosphere  temperature  and  composition.  Curves  are  drawn  for  blackbody 
temperatures  of  200,  220.  250  and  300  K.  Typical  radiation  temperatures 
between  5  and  25 yxm  lie  between  190  and  310  K  (E.  G.  Y/ark,  et  al. ,  1962). 

The  main  feature  of  importance  in  these  curves  is  the  decrease  in  vibrational 
temperature  with  increasing  wavelength  for  a  given  blackbody  temperature. 


During  the  day,  solar  infrared  flux  is  important  in  exciting  the  shorter 
wavelength  transitions.  Figure  III- 3  shows  solar  infrared  fluxes  at  the  top 
of  the  atmosphere  given  by  Johnson  (19645  and  the  resulting  vibrational 
temperatures  for  a  molecule  at  the  top  of  the  atmosphere.  Absorption  of 
radiation  by  molecules  between  the  lower  atmosphere  radiation  source  or  the 
sun  decrease?  the  effective  incident  radiation.  This  is  partially,  but  not 
completely,  offset  by  the  radiation  of  these  molecules. 


The  band  functions  are  used  to  determine  the  effects  of  absorption  between 
a  molecule  at  a  given  altitude  and  the  various  sources  of  radiation  affecting  a 
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given  band  in  the  following  way.  The  vertical  optical  paths  to  the  lower 
boundary  and  to  each  level  in  the  atmosphere  are  computed  and,  for  daytime 
conditions,  the  optical  path  to  the  assumed  top  of  the  atmosphere  in  the 


direction  of  the  sun.  The  effective  radiation  from  the  lower  boundary  is 
determined  by  multiplying  by  the  value  of  the  band  function  (  K^) 
corresponding  to  the  line  function  of  Equation  (3-43).  The  effective  radiation 
from  the  sun  is  determined  by  multiplying  by  the  band  function  M^(  KM 
corresponding  to  the  line  function  of  Equation  (3-42).  Finally,  the  effective 


radiation  from  all  atmospheric  levels  is  obtained  by  integrating  the  product 
of  emission  at  each  level  and  the  band  function  N^*(  K  ),  which  corresponds 


2  *  "o' 

to  the  line  function  of  Equation  (3-45).  Vibrational  populations  are  determined 
by  balancing  collisional  and  radiative  excitation  and  de-excitation  rates.  Total 
band  intensities  at  an  observation  point  inside  or  outside  the  atmosphere  are 
computed  by  determining  the  optical  path  from  each  level  to  the  observation 


point  and  integrating  the  product  of  local  emission  rate  and  the  band  function 
Kq)-  The  spectral  distribution  of  radiation  is  determined  by  finding  a 
mean  effectiv  .  temperature  and  line  of  sight  column  density  for  the  radiating 


molecules,  computing  intensity  for  each  line  in  the  band  with  the  line  function 


L(  frQ).  summing  the  line  emission  rates  and  scaling  each  line  emission  rate 


so  that  the  total  band  intensity  is  obtained. 


The  use  of  band  radiative  transfer  functions  derived  by  assuming  non¬ 
overlapping  lines  with  Doppler  shape  does  not  lead  to  a  completely  satisfactory 
integration  with  the  radiation  field  at  lower  altitudes.  Gille  and  Ellingson  (1968) 
have  computed  correction  factors  to  make  possible  the  inclusion  of  Doppler 
broadening  effects  in  the  computation  of  transmissivities  of  random  exponential 
bands  based  on  Lorentz  line  profiles.  A  rational  approximation  to  their  table 
of  correction  factors,  obtained  by  a  least  squares  gradient  search  on  the 
parameters  involved,  has  been  investigated.  Results  were  not  good,  partly 
because  of  the  inaccuracy  in  computing  transmissivity  differences  between 
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large  values  of  absorber  amount  and  partly  because  of  the  crudeness  of 
the.  rational  approximation,  which  in  some  cases  led  to  differences  with 
the  wrong  sign.  Another  problem  with  using  such  models  is  the  difficulty 
of  properly  accounting  for  irradiation  from  an  extended  source  (cf.  Kuhn 
and  London,  1969).  Line-by-line  integration  with  use  of  a  Voigt  line  profile 
as  done  by  Arnold,  et  al.  (1969)  is  feasible  when  vibrational  excitation  is 
known  beforehand  and  would  be  better  but  much  slower  than  the  present 
method  of  computing  limb  radiance  spectral  profiles. 

Table  III-l  lists  the  molecular  bands  included  in  the  program.  The 
band  strengths  listed  for  CO^,  H^O,  and  N^O  were  adopted  from  a 
compilation  by  C.  R.  Calfee  of  £SSA  which  was  provided  by  R.  A.  McClatchey 
(private  communication,  1970)  of  AFCRL.  These  total  band  strengths  agree 
with  those  obtained  by  summing  the  individual  line  strengths  used  by 
McClatchey,  et  al.  (1970).  The  ba-- ’  trength  of  the  fundamental  of  nitric 
oxide  is  an  average  of  the  values  o"  •  ed  by  Abels  and  Shaw  (1966),  Feinberg 
and  Camac  (1967),  Oppenheim,  et  <•  1967),  Schurin  and  Clough  (1963), 

and  Varanasi  and  Penner  (1967).  The  nitric  oxide  overtone  value  is  that 
obtained  by  Schurin  and  Ellis  (1966)  and  agrees  with  the  theoretical  calculation 
of  Michels,  et  al.  (1970). 

D.  Chemiluminescence 

It  is  well  known  that  some  reactions  such  as  that  of  atomic  hydrogen 
with  ozone  to  form  OH  lead  to  product  molecules  which  are  vibrationally 
excited  and  radiate  in  the  infrared  spectral  region.  Thus,  Corbin,  et  al. 

(1969)  assumed  that  nitric  oxide  produced  by  reaction  R42  was  produced 
with  an  average  of  one  vibrational  quantum.  Hushfar,  et.  al.  (1971)  have 
determined  the  initial  vibrational  populations  of  nitric  oxide  in  states  higher 
than  the  first  excited  level.  Their  results  indicate  that  at  most  one  molecule 
in  50  is  formed  vibrationally  excited  to  the  second  or  higher  vibrational  level. 


I 
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Table  HI-1 


Molecular  Band  Strengths  at  296K 


Vibrational 

Transition* 


Band  Center 

Upper 

Lower 

Band  Strength 

Molecule 

Wavelength,  jjjm. 

State 

State 

atm“^cm~2 

CO 

18.37 

11102 

10001 

0.00675 

U 

16.74 

11102 

02201 

0.  129 

16. 18 

10002 

01100 

3.56 

15.46 

11102 

10002 

0.  550 

14.98 

01100 

00000 

205.0 

14.  98 

02200 

01100 

16.  1 

14. 97 

03300 

02200 

0.948 

14.52 

11101 

10001 

0.369 

13.87 

10001 

01100 

4.  59 

13.48 

11101 

02200 

0.  196 

12.64 

11101 

10002 

0.0278 

10.41 

00010 

10001 

0.0122 

9.40 

00010 

10002 

0.0157 

5.  18 

11102 

00000 

0.0102 

4.82 

11101 

00000 

0.0553 

4.  26 

00010 

00000 

2380.0 

Hz° 

6.27 

010 

000 

290.0 

o3 

9.60 

001 

000 

320.0 

9.06 

100 

000 

8.26 

NO 

5.33 

1 

0 

112.0 

2.68 

2 

0 

1.95 

N.O 

17.26 

0200 

0110 

1.73 

o 

16.98 

0110 

0000 

30.4 

16.98 

0220 

0110 

1.73 

8.56 

0200 

0000 

11.  1 

7.78 

1000 

0000 

216.0 

7.74 

1110 

0110 

7.75 

4.58 

0001 

0000 

1710.0 

^Vibrational  states  are  designated  by  V,  V,  V,  ,  y  V7J?y,  or 

w^ere  V.  is  the  quantum  number  of  the  i^1  vibrational 
mode,  SL  is  the  vibrational  angular  momentum  quantum  number  and  f 
(if  greater  than  zero)  identifies  different  components  of  a  Fermi  resonance 
multiplet. 
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If  their  fit  of  the  distribution  to  a  relative  vibrational  temperature  of 
5000  K  holds  for  the  first  vibrational  level  also,  then  this  reaction  will  be  of 
negligible ‘importance  in  producing  infrared  radiation.  No  experimental  data 

i 

exists  on  the  vibrational  state  of  nitric  oxide  formed  by  the  more  important 
2 

reaction  of  N(  D)  with  molecular  oxygen,  and  nitric  oxide  chemiluminescence 
is  not  included  in  the  radiance  model.  Similarly,  F.  R.  Gilmore  (private 
communication  to  J.  W.  Carpenter,  1970}  has  suggested  that  ozone  may  be 
formed  in  highly  excited  vibrational  states  which  may  lead  to  radiation 
occurring  in  the  mode  at  longer  wavelengths  than  the  001  to  000  transition. 
Again,  there  has  been  no  experimental  determination  and  refinement. 
Gilmore's  estimates  woitld  require  more  information  on  both  the  initial  state 
of  the  product  ozone  molecule  and  the  relative  transition  probabilities  of  die 
competing  radiative  transitions  that  cascade  to  the  ground  vibrational  state. 
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SECTION  IV 


DESCRIPTION  OF  PROGRAMS 


This  section  describes  the  computer  programs  developed  to  implement 

:  i 

the  chemistry  discussed  in  Sfection  II  and  the  computation  of  high  altitude  limb 
viewing  infrared  radiances  employing  the  physical  processes  discussed  in 

i 

Section  in.  The  chemistry  program,  including  the  effects  of  vertical  eddy 


mixing,  and  molecular  diffusion,  is  completely  new.  The  programs  to  compute 


radiance  are  modifications  of  those  developed  by  Corbin,  et  al.  (1969).  An 


additional  program  has  been  written  to  selectively  plot  the  output  of  the  radiance 
program.  The  programs  are  described  below  and  are  listed  in  the  Appendix. 


A.  Chemistry  Program 


The  chemistry  program  consists  of  a  main  program,  SNAPS2, 
and  6  subroutines;  ROMTST,  TKSL,  SOLLY,  CONCEN,  COLUMN,  and 
SUNPHI.  The  main  program  f?riorms  the  chemical  computations,  including 
the  effects  of  vertical'  transport.  EOMTST  embodies  the  extrapolation  method 
of  Equation  (2-18).  TKEL  is  used  to  compute  temperatures  at  each  altitude. 
SOLLY  computes  the.  solar  zenith  angle  ;  nd  fraction  of  the  solar  disk  visible 
at  each  altitude.  COLUMN  integrates  the  number  densities  of  molecular 
nitrogen,  oxygen  and  ozone  above  a  given  altitude  in  the  direction  of  the  sun. 
CONCEN  is  called  by  COLUMN  to  obtain  number  densities  required  for  the 
integration.  SUNPHI  .computes  the  attenuation  of  solar  flux  by  absorption 
between  a  given,  altitude  and  the  assumed  top  of  the  atmosphere  in  the  solar 

*  i 

direction  and  returns  the  phbtodissociation  rates  required  for  the  chemistry. 
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The  necessary  input  data  to  the  chemistry  program  is  a  set  of 

14  data  cards  containing  identification  and  variable  parameters  and  a  set 

3 

of  603  cards  containing  concentrations  ofN,,  0_,  0(  P),  0_,  OH,  H,  H  O, 

1  *4  *  5  2  C 

H02,  H2Oz,  H2,  CH4,  C02,  CO,  0(  D),  N(  S),  NO,  N02»  N20,  N(  D)  and 
Ar  at  one  km  intervals  between  ground  level  and  200  km.  The  data  is  read 
in  the  following  order: 

1.  IDENT  60  characters  of  identifying  comments  which 

are  used  in  output  page  headings.  The  first 
character  should  be  a  blank. 

2.  DTHETA  Co-latitude,  degrees,  of  the  computation 

location.  A  value  of  zero  corresponds  to  the 
north  pole. 

3.  DPFI  East  longitude,  degrees,  of  computation 

location. 

4.  SOLDEC  Solar  declination  angle,  degrees,  used  to 

determine  time  of  year. 

2 

5.  GREF  Surface  gravitational  acceleration,  cm/ sec  , 

at  computation  location. 

6.  REFF  Effective  value  of  earth  radius,  km,  for  use 

in  computing  gravitational  accelerations  at 
points  above  ground  level. 

7.  EPS  Relative  error  permitted  for  each  time  step 

in  the  integration. 

8.  MLMAX  Maximum  number  of  times  integration  is  to  be 

attempted  for  each  time  step  to  obtain  relative 
error  EPS. 


IS- 


9.  NBASE  Index  of  base  altitude  for  computations, 
altitude  in  km  plus  one. 

10.  NTOP  Index  of  highest  altitude  in  computations, 

altitude  in  km  plus  one,  restricted  to  201  or 
less  in  current  program. 

11.  IHRMAX  Number  of  times  at  which  complete  output  is 

required. 


12.  MMAX 


Number  of  integration  steps  between  complete 
output  points. 


13.  DLTIME  Integration  interval,  seconds;  the  total  simulated 
time  for  a  given  run,  in  seconds,  is 
IHRMAX*MMAX*  DLTIME . 


14.  TIME 


15.  CONSP 


Initial  time,  seconds,  measured  from  noon; 
A.M.  values  are  negative,  P.M.  values  are 
positive. 

The  603  cards  containing  concentrations  and 
identifying  indices;  concentrations  are  read 
into  array  CONSP. 


Frogram  execution  begins  with  the  initialization  of  constants  used 
by  the  program.  The  values  of  eddy  diffusion  coefficients  uced  at  altitudes 
above  40  km  are  stored  in  array  EDCOEF.  The  peak  vaiue  set  by  the  DATA 
statement  occurs  at  104  km.  The  eddy  diffusion  coefficient  profile  is  lowered 
by  9  km  in  the  DO  loop  ending  with  statement  90  and  the  low  altitude  eddy 
diffusion  coefficients  in  array  EDCF  are  transferred  to  array  EDCDEF  in 
the  DO  loop  ending  with  statement  95.  This  permits  use  of  separate  sets  of 
high  and  low  altitude  diffusion  coefficients  by  changing  the  appropriate  data 
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statements,  and  easy  variation  of  the  altitude  of  maxium  turbulent  transport. 

The  data  described  above  is  then  read  in.  The  information  on  each  data  card 
is  written  to  printer  output  channel  6  as  it  is  read.  Columns  13  through  80 
of  the  serjnd  through  fourteen  data  cards  are  printed  in  addition  to  the  number 
placed  in  columns  2  through  12  so  that  comments  may  be  added.  The  first 
14  data  cards  are  printed  on  a  page  numbered  -1.  The  input  concentrations 
require  about  10  pages  to  print,  the  first  of  which  is  numbered  page  0.  Further 
execution  is  determined  by  the  input  data.  A  brief  flow  chart  is  given  on  the 
following  pages. 

The  outermost  control  DO  loop  ending  with  statement  450  is  executed 
IHRMAX  times.  Within  this  loop,  the  computer  clock  is  called  and  computer 
running  time  since  the  beginning  of  the  program  is  stored  in  variable  SECA. 
Variable  locations  used  in  the  DO  loop  ending  with  statement  110  are  set  to 
zero  and  this  loop  is  entered  and  executed  for  each  altitude  between  NBASE  and 
NTOP.  The  mean  molecular  mass  RMWT  at  each  altitude  is  computed  by 
the  statements  associated  with  DO  loop  102.  Temperatures  are  computed  by 
calls  to  function  TKEL.  The  gravitational  acceleration  at  each  altitude  is 
computed  and  local  mixing  scale  heights  are  computed.  Coefficients 

for  the  eddy  diffusion  contributions  to  q  and  terms  in  Equation  (2-10)  are 
determined  from  the  finite  differerce  approximation.  Equation  (2-8),  and  stored 
in  av.*ays  BTED,  QJEDL,  and  QJEDU.  Factors  needed  for  computation  of 
molecular  diffusion  contributions  are  computed  and  stored  in  arrays  AFAC, 
BFAC,  CFAC,  DFAC  and  RDTEMP.  The  molecular  concentrations  at  each 
altitude  above  the  base  altitude  are  scaled  to  be  consistent  with  the  finite 
difference  scheme  in  the  DO  loop  ending  with  statement  104.  An  effective 
molecular  scale  temperature  is  used  for  the  local  scale  height.  Before  exit 
from  DO  loop  110,  N0  and  O-  concentrations  and  sample  computed  values  are 

w  u 

written  to  the  printer  output  channel  for  altitudes  between  but  not  including  those 
associated  with  NBASE  and  NTOP. 
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The  DO  loop  ending  with  statement  400  is  entered  and  executed 
MMAX  times.  Each  execution  advances  the  integration  time  by  DLTIME 
seconds.  The  time  at  the  end  of  the  step  is  computed  in  floating  and  fixed 
point  form  If  the  value  of  NPRNT  permits,  a  new  page  is  started  and  the 
time  and  column  headings  are  printed.  Molecular  diffusion  coefficients  are 
computed  for  the  two  lowest  altitudes  in  the  statements  associated  with  DO 
loop  130. 

The  DO  loop  ending  with  statement  350  is  entered  and  executed  for 
each  altitude.  Molecular  diffusion  coefficients  are  computed  for  the  next 
higher  altitude,  and  the  molecular  diffusion  contributions  to  the  ^  and  q 
terms  are  computed  and  combined  with  the  eddy  diffusion  contributions  in 
the  DO  loop  ending  with  statement  140.  The  solar  zenith  angle  CHI  and  the 
fraction  ETA  of  the  sun's  disc  visible  are  computed  by  subroutine  SOLLY. 

If  the  sun  is  visible,  column  counts  of  N^»  G^,  O  and  O^  in  the  solar  direction 
are  computed  in  subroutine  COLUMN  and  photodissociation  rates  are  computed 
in  subroutine  SUNPHI  for  the  time  at  the  midpoint  of  the  time  interval  of 
integration.  Reaction  rate  coefficients  are  computed  for  the  temperature  at 
the  given  altitude. 

The  integration  for  a  global  time  increment  DLTIME  at  a  single  altitude 
is  performed  in  the  DO  loop  ending  with  statement  300.  The  operations  within 
this  loop  are  repeated  either  until  convergence  to  relative  accuracy  EPS  is 
achieved  or  until  MLMAX  repetitions  have  been  performed  with  no  convergence. 
Each  time  through  this  loop,  a  new  number  of  local  time  intervals,  MMF,  is 
computed  to  conform  to  the  extrapolation  procedure  of  subroutine  ROMTST. 

|  Values  of  the  initial  concentrations  at  the  integration  altitude  in  array  CONSP 

I  are  placed  in  array  DCX.  The  DO  loop  ending  with  statement  290  then  integrates 

[  over  the  MMF  local  time  steps.  The  values  in  array  DCX  are  copied  into  array 

|  TCX  and  the  time  incremer/t  DTIME  is  set  to  one-half  the  length  of  the  local 

j  increment.  An  index  ITER  is  set  to  one  and  appropriate  third  body  concentrations 
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are  computed.  The  algorithm  of  Equation  (2-10)  is  then  applied  for  each  species 
in  the  DO  loop  ending  with  statement  255,  the  results  being  placed  in  array  RCX. 
ITER  is  incremented  and  the  results  in  array  RCX  are  placed  in  array  TCX. 

The  time  increment  is  doubled  and  the  DO  loop  ending  with  statement  255  is  again 
executed,  using  the  intermediate  results  from  the  first  execution  to  compute 
the  aC  ,  ^3  and  q  terms.  The  results  of  this  second  cycle,  computed  in  array 
RCX,  are  transferred  to  array  DCX  and  the  statements  in  the  DO  loop  ending 
with  statement  290  repeated  until  the  global  time  interval  DL.TIME  has  been 
covered.  Subroutine  ROMTST  is  called  to  test  the  values  in  array  DCX  with 
previous  results  produced  in  the  DO  range  ending  with  statement  290.  The 
improved  concentration  values  are  returned  in  array  TCX.  If  convergence  is 
achieved,  control  is  transferred  to  statement  310.  Otherwise,  the  DO  loop  ending 
with  statement  300  is  repeated.  If  this  process  does  not  result  in  convergence 
in  MLMAX  tries,  the  endpoint  values  in  array  RCX  are  placed  in  array  TCX  in 
the  DO  loop  ending  with  statement  315.  The  endpoint  results  for  the  next  lower 
altitude,  which  have  been  stored  in  array  SCX,  are  written  into  array  CONSP 
and  the  contents  of  array  TCX  are  placed  in  array  SCX.  Production  rates  for 
OH  by  two  different  reactions  are  computed  and  stored  in  QCOHA  and  QCOHB. 
"Rectangular”  sums  of  total  column  counts  are  incremented  in  DZOHA  and 
DZOHB.  These  values  are  computed  to  provide  a  rough  check  on  the  hydrogen 
compound  chemistry  and  may  be  compared  with  airglow  observations.  If  the 
value  of  NPRNT  permits,  the  altitude,  QCOHA,  QCOHB,  DZOHA,  DZOHB  and 
the  concentrations  of  16  molecules  are  written  to  the  printer  output  channel. 

The  statements  in  the  DO  loop  ending  with  statement  350  are  then  repeated  until 
concentrations  at  all  altitudes  between  (but  not  including)  NBASE  and  NTOP 
have  been  integrated  over  the  global  time  step  DLTIME. 

When  all  altitudes  have  been  integrated,  a  summary  of  number  of 
iterations  is  sent  to  the  printer  output  channel  and  the  time  value  TIME  is 
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incremented  by  DLTIME. 

After  the  integration  haa  been  advanced  MMAX  times  with  the  increment 
DLTIME,  all  concentrations  are  sent  to  the  punch  output  channel  by  the  state* 
ments  in  DO  loop  420  and  to  the  printer  channel  by  the  statements  in  DO  loop 
440.  The  internal  computer  clock  is  again  accessed  and  the  starting,  end, 
and  elapsed  time  required  in  DO  loop  450  are  sent  to  the  printer  output  channel. 
Program  execution  ceases  after  this  outer  DO  loop  has  been  executed  IHRMAX 


times. 


The  scaling  of  concentrations  performed  in  the  DO  loop  ending  with 


statement  104  permits  use  of  output  data  from  runs  at  another  latitude  or  with 
a  different  temperature  profile  to  be  used  as  input  data.  While  there  may  be 
some  difficulty  in  reaching  convergence  at  altitudes  below  90  km  in  the  first 
time  step,  this  is  not  a  serious  problem.  When  using  a  global  time  step 
DLTIME  of  600  seconds  and  a  value  of  0.01  for  EPS,  computation  of  chemical 
changes  over  a  day  for  150  altitudes  requires  on  the  order  of  20  minutes  running 
time  on  the  CDC  6600.  This  is  longer  than  would  be  required  by  the  differencing 
methods  usually  employed  (e.g. ,  Shimazaka  and  Laird,  1970)  but  the  error 
control  compensates  for  this.  It  is  not  necessary  to  alter  time  increments  at 
dawn  or  sunset  since  the  convergence  algorithm  of  subroutine  ROMTST 
automatically  forces  the  use  of  smaller  time  subdivisions  of  the  global  time 
interval  DLTIME.  Further  study  is  required  before  it  will  be  possible  to 
state  optimum  values  of  DLTIME  and  MLMAX.  Since  as  many  as  20  diurnal 
cycles  may  be  required  to  reach  concentrations  which  repeat  from  one  simulated 
day  to  the  next,  a  means  of  reducing  running  time  is  required.  Possible  methods 
are  to  rewrite  the  program  and  subroutines  COLUMN,  CONCEN  and  TKEL  to 
use  a  larger  altitude  increment  or  to  generate  a  starting  atmosphere  with  the 
method  of  Colegrove,  et  al  (1966).  In  its  present  form,  the  chemistry  program 
SNAPS2  uses  fixed  concentrations  of  NO  ,  O^  ,  O  and  electrons  during  the  day 
and  assumes  no  ions  to  be  present  at  night.  The  full  ionic  chemistry  of  Table 
2-4  should  be  added  and  photoionization  rates  added  to  the  photodis60ciaiion 


rates  computed  by  subroutine  SUNPHI. 


The  names  of  the  real  arrays  and  their  uses  are: 


AFAC 

AMWT 

BFAC 

BTED 

BTEMD 

CARD 


stores  values  of  numerical  approximation  to 
_1  JT 
T  %z 

stores  values  of  molecular  weights,  gm/mole; 
values  are  set  by  data  statement 

stores  values  of  the  coefficient  of  molecular 
concentration  in  the  numerical  approximation 


stores  values  of  coefficients  used  in  computing 
contribution  to  terms  from  eddy  diffusion 

stores  values  of  contribution  to  term  from 
both  eddy  and  molecular  diffusion  for  each 
molecule  at  a  given  altitude  during  integration- 

used  as  temporary  storage  for  comments  in 
columns  13  through  80  of  first  14  data  cards 


CFAC 


stores  values  of  numerical  approximation  used 
1 

to  compute  — —  — —  for  diffusive  scale  heights 


CMD  stores  values  of  molecular  diffusion  coefficients 

for  each  molecular  species  at  a  given  altitude 
during  integration 

CMDL  stores  values  of  CMD  at  next  lower  altitude 

CMDU  stores  values  of  CMD  at  next  higher  altitude 
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COL 


CONSP 


DCX 


DFAC 


EDCF 


EDCOEF 


QJEDL 


QJEDU 


QJEMD 


RATE 


RCX 


stores  values  of  integrated  column  counts 


computed  by  subroutine  COLUMN 


stores  values  of  molecular  number  densities; 


initial  values  are  read  in  from  cards;  values  are 


altered  at  each  integration  step 


stores  initial  number  densities  for  one  altitude 


during  local  time  step  in  integration 


stores  values  of  g/RT  used  in  computing 
atmospheric  scale  heights  and  contribution  of 


molecular  diffusion 


stores  values  of  eddy  diffusion  coefficients  for 


altitudes  below  40  km,  set  by  DATA  statement 


stores  values  of  eddy  diffusion  coefficients. 


initially  contains  high  altitude  values  set  by 


DATA  statement 


stores  coefficients  of  eddy  diffusion  contribution  to 
q  terms  from  next  lower  altitude 


stores  coefficients  of  eddy  diffusion  contributions 
to  q  terms  from  next  higher  altitude 


stores  values  of  contribution  to  q  term  from  eddy 
and  molecular  diffusion  for  each  molecular  species 


at  a  given  altitude  during  integration 


stores  values  of  photodissociation  rates  computed 
by  subroutine  SUNPHI 


stores  intermediate  values  of  number  densities  for  a 


single  altitude  comp'  ced  during  local  time  step  in 
integration 


\  1 

‘i  « 


4  i 
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RDTEMP 


TCX 


stores  values  of  temperature  coefficients  used 
in  computing  molecular  diffusion  coefficients 

i 

stores  values  of  number  densities  computed  at 
end  of  global  time  step  for  next  lower  altitude; 
values  eventually  replace  corresponding  numbers 
in  CONSP 

stores  intermediate  values  of  number  densities 
used  in  computing  oi,  and  q  :  values  during 
integration  of  local  time  step  for  a  single  altitude 

I 

stores  values  of  temperatures  at  each  altitude, 
values  are  computed  by  function  TKEL 

stores  values  of  total  electron  densities  at  each 
altitude,  values  are  set  by  DATA  statement 

stores  format  statement  for  computer  running 
time  printout,  set  by  DATA  statement 

l 

stores  values  of  NO+  concentrations  at  each 

i 

altitude,  set  by  DATA  statement 

storea  values  of  0+  concentrations  at  each 
altitude,  set  by  DATA  statement 

stores  values  of  O*  concentrations  at  each 
altitude,  set  by  DATA  statement 

used  as  workspace  by  subroutine  ROMTST 


The  namf  i  of  integer  arrays  and  their  uses  are: 

IDENT  stores  60  character  comment  read  from  first 

data  card,  which  is  used  in  page  headings 


TEMPK 


TETOT 


TFMT 


TINO 


TIQ2 


IDENT 


.t-i.  . 


IDNT  stores  image  of  last  20  characters  of  first 

!  data  card 

■  I 

LSP  stores  hollerith  constants  used  in  printing 

column  headings 
1 

MLPRNT  stores  data  on  number  of  iterations  required 

for  convergence  of  integration  at  a  given  altitude, 
and ‘(largest)  index  of  species  which  did  not 

4 

converge  in  MLMAX  iterations 

The  names  of  the  real  variables  and  their  uses  are: 

AFACNA  temporary  storage  of  value  of  AFAC  for  altitude 

at  which  integration  is  being  performed 

BFACNA  temporary  storage  of  value  of  BFAC  for  altitude 
at  which  integration  is  being  performed 

BTEDNA  '  temporary  storage  of  value  of  BTED  for  altitude 
at  which  integration  is  being  performed 

i 

BTSQ  square  of^g  term  used  in  integration 

I 

BTT  value  (partial  or  complete)  of  ^  term  used  in 

integration  for  a1  single  molecular  species 

CFACNA  temporary  storage  of  value  of  CFAC  for  altitude 

at  which  integration  is  being  performed 

! 

CHI  solar  zenith  angle,  radians,  computed  by 

subroutine  SOLLY 

t 

CKISAV  temporary  storage  for  value  of  CHI 

CMDSNP  •  temporary  storage  for  value  of  molecular  diffusion 

coefficient  for  single  molecular  species 
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CON 


CONFAC 


CONSPN 


Cl 


C2 

C3 

DCT 


DDDZ 


DELZ 


DFACNA 


DHDZ 

DKDZ 

DLT 


temporary  storage  for  factor  used  in  scaling 
molecular  concentrations  at  altitude  NTOP  to 
to  fit  scale  height  defined  by  molecular  scale 
temperature 

temporary  storage  for  factor  used  to  scale 
molecular  concentrations  to  fit  local  scale 
heights  defined  by  molecular  scale  temperature 

temporary  storage  for  concentration  of  a  single 
molecular  species  at  a  given  altitude 

temporary  storage  for  molecular  nitrogen 
concentration 

temporary  storage  for  molecular  oxygen 
concentration 

temporary  storage  for  atomic  oxygen 
concentration 

temporary  storage  for  single  value  from 
array  DCX 

temporary  storage  for  value  of  — 1 —  for  single 

t>  z 

molecular  species 
altitude  increment,  cm 

temporary  storage  for  value  of  DFAC  at  altitude 
for  which  int  egration  is  being  performed 

temporary  storage  for  value  of 

temporary  storage  for  value  of 

02 

value  of  £  used  in  integration  for  a  single 
molecular  species 


1 
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DLTIME 

DPHI 


DTDZ 

DTHETA 


DTIME 

DZOHA 

DZOHB 

ED 


EPS 

ETA 

FAC 


FACB 

FACDMB 
FACM 
FACRT 
FAC  2  AX 


time  increment  for  integration  step, read  in 

as  data 

east  longitude  of  point  of  integration,  read  in 
as  data;  not  needed  in  present  program  which 
assumes  times  are  true  solar  times 

JT 

temporary  storage  for  value  of  -r-* 

07, 

colatitude  of  point  at  which  integration  is 
performed,  read  in  as  data 

storage  for  value  of  locai  time  step  in  integration 

values  of  integrated  production  rates  of  OH 

temporary  storage  for  value  of  eddy  diffusion 
coefficient  at  a  given  altitude 

value  of  relative  error  allowed  in  integration 

fraction  of  solar  disc  visible  at  given  altitude 

temporary  storage  for  intermediate  results, 
use  varies  in  different  parts  of  program 

temporary  storage  for  either  yit  or  JSt  in 
Equations  (2-10) 

temporary  storage  for  Sye  in  Equation  (2-10a) 

temporary  storage  for  value  of  1  -exp(-  £  t) 

2 

temporary  storage  for  4-eCq/jB 
temporary  storage  for  2ocxq 


FAC4AQ 


temporary  storage  for  4«<q 


G 

GL 

GREF 

GU 

QCOHA 

QCOHB 

QJEDLN 

QJEDUN 

Q.TT 

RDELZ 

RE 

REFF 

RFAC 
RI1, . . . 
RMWT 


value  of  gravitational  acceleration  at 
altitude  NALT 

value  of  gravitational  acceleration  at 
altitude  NALT  -  1 

surface  value  of  gravitational  acceleration, 
read  in  as  data 

value  of  gravitational  acceleration  at 
altitude  NALT  +  I 

values  of  OH  production  rates  at  altitude 
of  integration 

temporary  storage  for  value  of  QJEDL  at  altitude 
of  integration 

temporary  storage  for  value  of  QJEDU  at 
altitude  of  integration 

temporary  storage  for  partial  or  complete 
value  of  q  term  used  in  integration  of  a  single 
molecular  species 

reciprocal  of  twice  the  altitude  inci  'ir.ent  DELZ 

effective  radius  of  earth,  km,  used  in  geometrical 
computations 

effective  radius  of  earth,  km,  used  in  gravitational 
calculations,  read  in  as  data 

temporary  storage  of  intermediate  results 

storage  for  values  of  ionic  reaction  rates 

mear  -  .olecular  mass  at  altitude  NALT 
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RMWTL 

RMWTU 

ROBS 

RTIME 

RL  •  •  • 

R2DELZ 

SECA 
SECB 
SECC 

SFAC 

SH 

SHL 

SHU 

SU,  ... 

SOLDEC 

SUM 

SUML 

SUMU 

SI,  ... 

TCH4,  .. 


mean  molecular  mass  at  altitude  NALT  *  1 

mean  molecular  mass  at  altitude  NAi_,T  +  1 

altitude  of  given  point,  km 

variable  use  1  in  printing  time  values 

values  of  forward  neutral  reaction  rates 

square  of  reciprocal  of  altitude  increment,  DELZ 


storage  for  value  of  computation  of  computer 
runni  ng  time 

temporary  storage  of  intermediate  results 
value  of  mixing  scale  height  at  altitude  NALT 
value  of  mixing  scale  height  at  altitude  NALT  -  1 
value  of  mixing  scale  height  at  altitude  NALT  +  1 
values  of  reverse  ionic  reaction  rates 
solar  declination  angle,  read  in  as  data 


temporary  storage  used  in  forming  mean 
molecular  masses 


values  of  reverse  neutral  reaction  rate 
coefficients 

values  of  individual  molecular  concentrations  in 
array  TCX,  names  individually  for  convenience 
in  forming  cC ,  and  q  coefficients;  TCH4 
corresponds  to  CH.,  TCO  corresponds  to  CO, 

®r 

etc. 
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TEMP 

TEMPI 

TEMPL 

TEMPU 

TEMP12 

TEMP32 

TEMP52 

TFAC 

TIME 

XM4,  .  . . 

XXX12 

The  names 

I 

IF 


value  of  temperature  at  altitude  NALT 

reciprocal  of  temperature  1 

value  of  temperature  at  altitude  NALT  -  1 

value  of  temperature  at  altitude  NALT+  1 

value  of  minus  one-half  power  of  temperature 

value  of  minus  three-halves  power  of  temperature 

value  of  minus  five -halves  power  of  temperature  ■ 

temporary  storage  of.intermediate  results 

time,  expressed  in  seconds  from  noon,  read  in 
as  data;  altered  in  program  to  be  the  time  at 
midpoint  of  a  global  integration  step  so  th~t 
approximate  mean  values  of  photodissociation 
rates  will  be  computed  by  subroutine  SUNPHI 

third  body  concentrations  for  appropriate 
three-body  reactions 

storage  for  factor  used  in  computing  q  term 
for  C and ^  term  for  CO 

of  integer  variables  and  their  uses  are: 

I  ’ 

index  i"  DO  loop  95,  auxiliary  index  in  DO 
loop  110,  and  index  in  implied  input/output  : 

DO  loops 

final  parameter  in  innermost  implied  DO  loop 
used  in  printing  output  concentrations  within 
DO  loop  440 


IHR 

IHRMAX 

l 

IP 

) 

IS  ’ 

ITJER 


DC 


J. 


KCH4,  ... 


KTEST 


index  in  DO  loop  450 

final  parameter  in  DO  loop  450,  input  as 
data 

index  in  print  control  DO  loop  440 

initial  parameter  in  implied  DO  loop,  see  IF 

control  index  in  integration,  used  between 
statements  200  and  260  to  select  proper  values 

r 

for  array  TCX 

index  in  DO  loop  435,  incremented  for  each 
page  of  output 

index  in  implied  DO  loops  used  in  reading  data 
into  array  CONSP  within  DO  loop  100 

index  variables  used  to  select  appropriate 
values  of  quantities  related  to  given  species; 
most  uses  replaced  by  use  of  TCH4,  etc, 

parameter  returned  by  subroutine  ROMTST, 
zero  if  integration  has  converged,  non-zero 
otherwise 


LALT,  LCHI, 
LCOLN2, 
LCOL02, 
LCOL03, 
LDZOHA, 
LDZOHB, 

L.J02,  LJ03, 
LML,  LQCOHA, 
DQCOHB, 
LTEME 


store  Hollerith,  constants  used  in  printing  column 

i 

headings  and  are  set  by  DATA  statements 


1 

1 

Si  < 

1 

1 


1 

1 

5 
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MLMAX 


index  of  DO  loop  400,  incremented  for  each 
global  time  step 

index  of  DO  loop  300,  incremented  for  each 
iteration  of  integration  over  global  time  step 
at  a  given  altitude 

final  parameter  of  DO  loop  300,  read  in  as 


MLSAVE 


MMAX 


temporary  storage  for  value  of  ML,  used  after 
exit  from  DO  loop  300  to  print  number  of 
iterations  required  for  convergence 

final  parameter  of  DO  loop  400,  read  in  as 


MMF 


MMS 


NALT 


NAM  IN 


NBASE 


NMXN 


final  parameter  of  DO  loop  290,  computed  from 
ML  and  is  number  of  local  subdivisions  of  global 
time  step  DLTIME  for  the  MLth  iteration  of 
integration 

index  for  DO  loop  290,  see  MMF 

auxiliary  index,  usually  set  to  NALT  -  1 

index  for  DO  loops  90,  100,  110,  350,  420  and 
430;  denotes  a  given  altitude 

value  of  NALT  -  1  in  DO  loop  350 

value  of  NALT  -  2  in  DO  loop  110 

index  of  lowest  altitude  for  integration,  read 
in  as  data 

final  parameter  of  DO  loop  430 

stores  input  value  of  NBASE,  which  is  changed 
by  program 
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T 


NPAGE 


NPLU 


NPRNT 


NSPM 


NTIME 


NTOP 


NXX 


stores  output  page  number 


value  of  NALT+  1  in  DO  loop  350 


variable  available  to  control  output  of  concentrations 
at  integration  steps  intermediate  to  fully  documented 
output  at  end  of  DO  loop  450,  set  to  1  at  present,  so 
values  are  output  for  each  interval  DLTIME 


initial  parameter  of  DO  loop  430 


index  in  DO  loops  102,  104,  130,  140,  180,  190, 


255,  258,  280,  315,  and  320;  a  value  of  NSP 


selects  a  value  of  a  variable  corresponding  to  a 
single  molecular  species 


number  of  molecular  species  included  in  program, 
20  at  present,  used  as  final  parameter  of  DO  loops 


with  NSP  as  index 


final  parameter  of  DO  loop  420,  in  which  output 
cards  are  punched 


used  in  printing  time  values 


index  of  highest  altitude  used  in  integration, 
input  as  data 


channel  number  for  printing  computer  running  time 
information,  set  to  6  here 


temporary  storage  for  altitude  index  read  from 
cards  containing  concentration  data 


storage  space  used  for  sequence  numbers  of  cards 
read  in  or  punched  out  at  a  given  aldtude 


B.  Radiance  Program 


The  high  altitude  radiance  programs  of  Corbin,  et  al.  (1969), 
BCKGND  and  SPCTRA,  have  been  modified  by  the  addition  of  six  carbon  dioxide 


i 

! 

! 


I 


bands  between  5  and  25 yMjm,  and  improved  band  modeling  for  the  water  vapor 
6.3/<un  band  and  the  ozone  9.  1  and  9.6yLtm  bands-  The  basic  structure  is 
unaltered  and  the  division  into  two  programs  is  retained.  The  organization 
of  the  CDC  6600  program  control  system  makes  this  feasible  and  the  overall 
saving  in  central  memory  storage  space  is  greater  than  would  be  possible  by 
using  overlays.  The  output  from  the  three  lower  altitude  programs  of  Corbin, 
et  al  (1969),  which  are  based  on  Anding*s  (1968)  program,  is  required  only  to 
produce  values  of  vertically  directed  radiances  at  the  band  centers.  The 
programs  INTERP  and  CATM  were  run  to  produce  on  cards  a  complete  set  of 
input  data  to  RAD  for  each  of  the  47  lower  atmosphere  models.  RAD  was  modified 
slightly  to  permit  reading  in  the  band  data  cards  and  to  either  write  on  tape  or 
punch  the  band  data  cards  with  appropriate  values  of  upwelling  flux,  TFLUX. 

The  mid -latitude  winter  atmosphere  model  was  used  for  the  present  computations. 

Program  BCKGND  consists  of  the  main  program  BCKGND  and 
subroutines  BANRAD,  INTERP,  SIMP,  XFER,  TAUMAX,  SJMAX,  SDM2,  SMI, 
SLM1,  QJPART,  and  EVAL.  Program  SPCTRA  contains  the  main  program 
SPCTRA  and  subroutines  ROTATE,  SFUN,  QINT,  FILTER,  WATSTR,  ROTSTR, 
ROTPOS,  and  SIMP.  Overall  operation  is  discussed  below. 


I?> 


mXr 
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The  main  program  BCKGND  serves  primarily  to  control  the  operation 
of  subroutines  which  compute  vibrational  populations  and  radiative  rates. 
BCKGND  reads  and  checks  the  input  data,  and  transfers  control  to  subroutine 
BANRAD  which  computes  vibrational  populations,  local  optical  thicknesses 
and  volume  radiances  for  all  incl  uded  bands  of  a  single  molecular  species. 
When  control  is  returned,  BCKGND  computes  and  outputs  integrated  band 
radiances  for  limb  viewing  (exoatmospheric)  and  seven  endoatmospheric 
viewing  angles  at  selected  tangent  heights  or  endoatmospheric  altitudes. 

Data  is  read  in  on  two  different  channels,  one  being  used  for  number 
densities  and  collisional  excitation  rates,  the  other  for  band  data  which  is 
needed  again  by  program  SPCTRA.  The  storage  space  used  for  methane  data 
in  the  older  version  is  now  used  for  atomic  oxygen  concentrations.  After  input 


data  is  checked,  control  is  transferred  to  subroutine  BANRAD. 

BANRAD  computes  the  population  of  vibrational  levels  of  infrared 
active  species,  including  the  effects  of  collisional  exciation;  absorption  of 
radiation  from  the  lower  atmosphere,  sun  and  molecules  radiating  at  high 
altitude.  The  values  returned  to  the  main  program  are  volume  radiance  and 
unit  optical  thickness  at  each  altitude  for  each  band  included  in  the  computation. 

The  first  step  in  the  computation  is  determination  of  the  thermal 
equilibrium  populations  of  vibrational  levels.  This  is  performed  in  the  DO 
range  ending  with  statement  130.  Each  passage  through  this  DO  range  computes 
populations  for  one  altitude.  Within  this  DO  range  are  two  others.  The  DO 
range  ending  with  statement  110  computes  the  vibrational  population  of  each  level 
relative  to  the  ground  vibrational  level,  using  the  energies  (wavenumber)  stored 
in  WAVE  and  the  statistical  weights  stored  in  the  ten’s  places  in  LSC.  The  relative 
populations  are  stored  temporarily  in  real  array  VIBPOP.  The  sum  of  relative 
populations  is  stored  in  SUM.  The  final  ground  state  population  is  then  computed 
by  dividing  total  concentration  by  SUM,  and  stored  in  real  array  GNDSTE.  The 
DO  range  ending  with  statement  120  then  computes  the  vibrational  population  of 
each  level  by  multiplying  total  concentration  by  relative  configuration  divided  by 
SUM. 

The  DO  range  ending  with  statement  135  selects  the  vibrational  level 
which  is  coupled  with  Hie  first  excited  nitrogen  vibrational  level. 

The  DO  range  ending  with  statement  160  computes  the  rate  at  which  a 
single  molecule  in  a  given  level  is  excited  or  de-excited  by  collisions  at  each 
altitude. 

The  ovt  er  DO  range  ending  with  statement  180  computes  wavelengths 
(BNDLAM),  Einstein  A  coefficients  (EA),  G  factors  for  absorption  of  radiation 
from  the  sun  (GDOWN)  and  from  the  earth  and  lower  atmosphere  (GUP)  and  the 
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G  factor  used  in  computing  radiative  transfer  within  the  upper  atmosphere 
(GFAC).  The  inner  DO  range  ending  with  statement  180  computes  the  unit 
optical  thicknesses  for  each  band  at  each  altitude  using  function  TAUMAX, 
and  stores  the  result  in  array  BNDTAU.  In  addition,  it  computes  the  volume 
radiances  and  stores  them  in  array  BNDRAD.  Further  computation  is  then 
performed  iteratively  in  the  DO  range  ending  with  statement  400. 

The  DO  range  ending  with  statement  190  computes  various  integrated 
optical  thicknesses  for  each  band.  The  DO  range  ending  with  statement  182 
computes  the  total  optical  thickness  between  each  two  altitudes  in  the  vertical 
direction  and  stores  the  values  in  array  TTH.  The  DO  range  ending  with 
statement  184  (1)  computes  the  total  optical  thickness  between  the  base  altitude 
and  each  other  altitude,  storing  the  result  in  array  TTTT?  and  (2)  computes 
the  total  optical  thickness  between  the  top  altitude  and  each  other  altitude, 
storing  the  results  in  array  TTDOWN.  If  the  computation  is  not  for  daytime 
conditions,  the  DO  range  ending  with  statement  186  sets  the  values  in  TTDOWN 
equal  to  0.0,  for  efficiency  in  the  later  computations. 

Computation  of  vibrational  populations  is  done  separately  for  each 
species.  Control  is  transferred  to  statement  220  for  CO^,  statement  230 
for  H^O,  statement  240  for  NO,  statement  250  for  N^O  or  statement  230 
for  O^*  The  DO  range  following  e^ch  of  these  statements  is  used  to  compute 
vibrational  populations.  After  each  set  of  vibrational  populations  is  completed, 
new  values  of  volume  band  radiances  are  computed  in  the  nested  DO  ranges 
ending  with  statement  360.  Following  three  iterations,  vibrational  populations 
are  written  to  an  output  file  for  use  by  program  SPCTRA  and  control  is  returned 
to  the  main  program. 

BANRAD  returns  volume  emission  radiance  values  in  array  RAD  and 
local  optical  thickness  values  in  array  TAU.  For  each  band,  line  of  sight 
band  radiances  are  computed  in  the  DO  range  ending  with  statement  840.  Then 


trapezoidal  rule  integration  is  performed  to  determine  the  radiance  arising 
along  each  0|f  the  eight  lines  of  sight,  including  the  modification  necessary 

i  . 

forthe  effects  of  optical  thickness  along  the  line  of  sight.  The  DO  range  ending 
with  statement  830  controls  the  integration  for  each  tangent  height  and  endo- 
ktmospheric  viewing  level.  Radiance  values  are  computed  at  5  km  intervals 
of  tangent  height  up  to- 150  km  and  at  25  km  intervals  above  this  level.  The 
integration  of  radiance  values  is  performed  for  each  endoatmq spheric  viewing 
level  in  the  DO  range  ending  with  statement  810.  The  integration  begins  ai 
the  viewing  altitude.  At  each  point  on  the  line  of  sight,  the  contribution  of  the 
local  volume  emission  rate  to  that  received  at  the  viewing  point  is  computed 
by  subroutine  SMI  which  takes  into  account  the  effects  of  optical  thickness  between 
the  emitting  point  and  the  viewing  point.  The  equivalent  integration  for  the  limb 
viewing  case  is  performed  in  the  DOranges  ending  with  statements  815  and  818. 
Upon  completion  of  the  integration  for  each  tangent  height  and  endoatmospheric 
viewing  level,  bfind  radiances  for  that  level  are  both  printed  and  written  to  the 
output  file  used  by  BANRAD,  for  input  to  SPCTRA. 

'  Program  SPCTRA  computes  the  spectral  radiances  between  5  and  25y£<cm 

corresponding  to  the  band  radiances  computed  by  ECKGND.  The  input  data  to 

SPCTRA  are  the  molecular  vibrational  level  and  band  data  also  used  by  BCKGND, 

the  vibrational  population  and  band  radiance  data  generated  by  BCKGND,  and 

individual  line  data  for  the  water  vapor  rotational  lines  and  band  and  the 

ozone  V  and  V  bands.  The  later  data  was  obtained  from  R.  A.  McClatchey 
3  1  , 

through  the  contract  monitor. 

When  program  execution  begins,  values  are  assigned  to  the  integer 
i  variables  NOUT,  NLN  and  NAL.T  and  the  subroutine  ROTATE,  which  reads  Ln 

i 

data  for  the  individual  lines,  is  called.  The  outer  DO  range  ending  with  statement 
500  is  entered.  The  DO  parameters  are  fixed  at  present  to  produce  two  sets  of 
spectra,  using  night  and  day  input  from  program  BCKGND.  The  program  may  be 
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made  more  flexible  by  reading  in  DO  parameters  prior  to  execution  of  the  DO 
statement.  The  variable  BLAM  is  set  to  5.0  and  the  outer  DO  range  ending  with 
statement  105  is  entered.  The  effect  of  the  first  two  statements  is  to  set  up  a 
table  of  wavelengths  at  0.  1  juxn.  intervals  between  5  and  25  jjjm.  in  real  array 
AVLAM.  These  are  the  wavelengths  at  which  spectral  radiances  are  printed 
out.  Again,  more  flexibility  may  be  achieved  by  making  the  intervals  and 
spacing  dependent  on  values  read  in.  If  this  is  done,  corresponding  changes 
must  be  made  in  subroutine  FILTER.  The  inner  DO  range  ending  with  statement 
105  sets  the  values  of  output  spectral  radiances  equal  to  zero.  Temperatures 
are  set  in  the  DO  range  ending  with  statement  110. 

Next,  the  number  of  molecular  species  is  read  in  and  stored  in  integer 
variable  NSPTOT.  The  outer  DO  range  ending  with  statement  400  is  entered 
and  executed  NSPTOT  times.  Within  this  DO  range,  molecular  vibrational  level 
and  band  data  and  vibrational  population  data  are  read  in  and  mean  temperatures 
and  total  number  of  molecules  per  unit  area  in  the  vertical  and  horizontal 
directions  are  computed  for  each  vibrational  stc.te.  Total  band  radiances 
corresponding  to  limb  viewing  and  vertical  and  horizontal  endoatmospheric 
viewing  are  read  in.  Relative  line  radiances  are  computed  and  normalized  to 
the  total  band  radiances.  The  spectra  corresponding  to  the  normalized  line 
radiances  are  computed  and  stored  in  arrays  SPECA,  SPECB  and  SPECG.  Details 
of  thes  computations  are  given  below. 

The  code  number  of  the  species  ISP  is  read  in.  The  number  of  vibrational 
levels  NLEVEL,  the  number  of  bands  NBAND  and  the  molecular  weight  RM  are 
read  in.  In  the  DO  range  ending  with  statement  210,  the  level  description  codes 
LSC  and  LCC,  the  energy  and  the  mean  rotational  constant  for  each  vibrational 
level  are  read  in.  In  the  DO  range  ending  with  statement  220  the  band  description 
codes  LBC  and  LBU  and  the  band  strength  for  each  band  are  read  in.  The  number 
of  altitudes  NALT  for  which  vibrational  populations  are  computed  for  the  species 
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and  the  number  of  altitudes  NRLEV  for  which  band  radiances  are  computed 
by  BCKGND  are  read  in.  Vibrational  level  populations  are  read  in  and  stored 
in  array  CON. 

The  mean  temperatures  and  total  number  of  molecules  per  square 
centimeter  in  the  vertical  and  horizontal  directions  at  5  km  intervals  between 
60  and  115  km  are  determined  in  the  DO  range  ending  with  statement  129.  The 
product  of  temperature  and  vibrational  level  population  is  computed  in  the  DO 
range  ending  with  statement  112.  The  integrals  of  total  number  of  molecules 
and  the  product  of  temperature  and  number  density  in  the  vertical  direction  are 
computed  in  the  DO  range  ending  with  statement  120.  Simpson's  rule  integration 
is  used.  The  total  number  of  molecules  in  a  given  vibrational  level  above  a  given 
altitude  is  stored  in  array  COLCNT.  The  corresponding  mean  temperatures, 
obtained  by  dividing  the  integral  of  the  product  of  temperature  and  number 
density  by  the  intergral  of  number  density,  are  stored  in  array  AVTEMP. 

After  exit  from  DO  range,  the  column  counts  and  temperatures  are  printed  out. 

The  DO  range  ending  with  statement  128  performs  a  ?  'ar  computation 
for  column  counts  and  mean  temperatures  in  the  horizontal  direction.  The 
integrations  required  are  performed  in  function  QINT-  Upon  exit  from  this  DO 
range,  column  counts  stored  in  array  HORCNT  and  mean  temperatures  stored 
in  array  BVTEMP  are  printed  out. 

The  inner  DO  range  ending  with  statement  400  is  now  entered  and  executed 
for  each  band.  This  DO  range  controls  the  computation  of  the  spectra. 

The  altitude  HT,  the  limb  band  radiance  RADC,  and  endoatmospheric 
radiances  RADA  (upward  direction)  and  RADB  (horizontal  direction)  are  read 
in  the  DO  range  ending  with  statement  230.  Band  information  is  computed 
from  band  and  level  constants.  A  test  is  made  on  the  species  code  ISP  to 
determine  whether  the  molecule  is  linear  or  nonlinear. 


Following  the  test  of  ISP,  if  the  molecule  is  linear,  line  positions  are 

computed  in  subroutine  ROTPOS.  The  DO  range  ending  with  statement  250 

! 

is  entered  and  executed  for  each  of  the  twelve  altitudes  for  which  spectra  are  . 
computed.  Line  radiances  are  computed  in  subroutine  ROTSTR  for  Vertical 
endoatmo spheric  viewing.  Successive  calls  to  subroutine  FILTER  compute 

1  i 

the  spectral  contributions  from  lines  in  the  P,  Q,  and  R' branches.  The 
sequence  of  calls  to  ROTSTR  and  FILTER  is  repeated  for  endoatmo  spheric 
viewing  in  the  horizontal  direction  and  for  limb  viewing.  Upon  exit  from  this 

t  ^  i 

DO  range,  control  is  transferred  to  statement  28(L 

If  the  test  on  ISP  determines  that  the  molecule  is  non-linear  (water 
vapor  or  ozone)  control  is  transferred  tp  statement  260  and  the  DO  range 

!  *  i 

ending  with  statement  270  is  executed  for  each  of  the  altitudes  for  which  spectra 
are  computed.  For  water  vapor,  spectra  are'  generated  from  the  pure  rotational 
line  radiances  by  calles  to  subroutine  WATSTR  and  FILTER  and  for  the  6.  3^itm 
band  by  calls  to  WATSTR  and  FILTER,  for  each  of  the  three  viewing  conditions. 
For  ozone,  spectra  in  the  9.  1  and. 9- 6  ^un  bands  are  also  generated  by  calls 
to  WATSTR  and  FILTER. .  ! 

i  ! 

Following  the  final  execution  of  the  DO  ranges  ending  with  statement 
400,  the  spectra  are  printed  out  in  the  order  endoatmospheric  viewing  in  the 
upward  direction,  endoatmospheric  viewing  in  the  horizontal  direction  and 

i 

limb  viewing.  Each  of  these  sets  is  printed  on  four  pages,  successive  pages 
covering  5-10ytun,  10-15/&m,  15-ZO^itan  and  ,20-25//4tm.  The  first  column 
on  each  page  contains  wavelengths  and  successive  columns  contain  radiance 
values  at  altitudes  or  tangent  heights  60,  65,  70  ,  ...  km.  The  data  in  : 

*  i 

arrays  AVLAM,  SPECA,  SPECB  an$3  SPECC  are  written  to  a  binary  file, 
either  for  immediate  tise  by  the  plotting  program  or  for  transfer  to  magnetic 

i 

tape  for  later  use. 
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C.  Plotting  Program  1 

The  plotting  program,  TRYCRT,  is  designed  to  read  the  binary 
output  tape  produced  by  the  spectral  radiance  program,  SPCTRA.  It  consists 
of  the  main  program,  TRYCRT,  and  six  subroutines,  I.INSCA,  LOGSCA,  LENAXS, 

i 

LOGAXS,  GIUD  and  NULINE.  In  addition,  it  uses  the  AFCRL/CDC  6600  system 

I 

.  subroutines  which  are  needed  to  produce  microfilm  plots.  The  version  listed 
in  the  Appendix  is  se.t  up  to  plot  three  sets  of  radiances,  assuming  that  the 
SPCTRA  output  tape  contains  two  sets  of  data,  with  vertical  and  horizontal 
endoatmospheric  radiances  and  limb  radiance  data  in  each  set.  This  program 
further  assumes  that  there  are  12.  different  altitude  profiles  in  each  data  set* 
Graphic  output  is  written  on  tape  39  and  microfilm  is  produced  from  this.  The 
microfilm  must  then  be  reproduced  on  paper. 

The  program  places  the  date  of'  the;run  in  the  first  work  of  array  IDENT, 
which  is  plotted  with  each  graph.  The  other  2  words  of  the  array  contain  the 
user  name  and  project  number,  stored  by  a  data  statement.  The  x-axis  length, 
SX,  is  set  to  9  inches  and  the  y-axis  length,  SY,  is  set  to  14  inches.  The  CDC 
6600  subroutines  are  initialized  by  calls  to  subroutines  CRTPLT  PLOT. 


f 


The  number  of  graphs  to  be  plotted,  NCASES,  number  of  characters^  NXL,  in 
x-axis  caption,  k-axis  caption,  XBCD,  number  of  characters,  NYL,  in  y-axis 
caption,  and  y-axis  caption,  YBCD,  are  read  in  from  cards  and  printed  on  the 
line  printer.  'The  DO  loop  ending  with  statement  400  is  then  executed  for  each  ! 
graph  to  be  produced.  Minimum  wavelength,  WVMIN,  and  maximum  wavelength, 
WVMAX,  are  read.  Plotting  constants  INF,  NX  and  NY  are  computed.  Wave¬ 
lengths  are  read  from  the  SPCTRA  output  tape  (channel  2)  and  stored  in  array 
AVLAM.  An  appropriate  «ubset  is  then  stored  in  array  X.  These  are  then: 
scaled  by  subroutine  LIN$CA.  Radiance  data  are  then  read  from  channel  2  and 
stored  in  arras  SPCTRA.  An  appropriate  subset  is  selected  and  stored  in  array 
Y.  Values  in  Y  are  scaled,  for  plotting  by  a  call  to  subroutine  LOGSCA.  Line 
intensity  is  set  to  16  by  a  call  to  subroutine  STBEAM.  Axes  Are  drawn  on  all 
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four  sides  of  the  graph  by  calls  to  subroutines  LINAXS  and  LOGAXS.  Calls 
to  subroutine  GRID  place  vertical  and  horizontal  grid  lines  on  the  graph. 

Radiance  values  are  further  scaled  in  the  DO  loop  ending  with  statement  150. 
Radiance  values  are  placed  on  the  graph  by  calls  to  subroutine  NULINE  in 
the  loop  ending  with  statement  200.  Line  intensities  for  each  curve  produced 
in  this  loop  are  set  by  calls  to  ST3EAM  and  use  values  stored  in  array  JBEAM 
through  a  data  statement.  Alter  each  set  of  curves  on  a  graph  is  produced, 
a  calle  to  subroutine  FRAME  piepares  for  a  new  graph.  The  DO  loop  ending 
with  statement  250  is  repeated  three  times,  the  first  and  second  cycles  plotting 
vertical  and  horizontal  endoatmospheric  radiances  while  the  third  plots  limb 
viewing  radiances.  When  only  limb  viewing  radiances  are  desired,  statement 
250  should  be  placed  immediately  after  the  statement  "read  (2)  SPCTRA" 
on  line  TCRT  540. 

Subroutines  LINSCA  and  LOGSCA  were  written  to  prov.de  automatic 
scaling  of  data  points  in  a  manner  which  would  utilize  maximum  space  on  the 
graph.  In  addition  to  scaling  data  points,  these  subroutines  generate  values 
KONX,  DX,  MINY  and  MAXY  which  provide  data  which  is  used  in  the  axis 
drawing  subroutines  LINAXS  and  LOGAXS  and  the  grid  drawing  subroutine 
GRID.  The  subroutine  NULINE  was  written  to  replace  the  standard  AFCRL/CDC 
6600  subroutine  LINE  which  draws  straight  lines  between  plots.  NULINE 
permits  points  to  be  plotted  a  limited  distance  outside  the  graph  boundaries 
if  the  data  for  a  graph  covers  more  than  9  orders  of  magnitude. 
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SECTION  V 


RESULTS  AND  DISCUSSION 


Representative  limb  viewing  radiance  profiles  are  shown  in 

Figures  V-l  through  V-5.  These  profiles  result  from  the  use  of  the 

nitric  oxide,  water  vapor,  ozone  and  carbon  dioxide  concentrations  of 

Figures  II- 1  to  Q-4  and  an  assumed  nitrous  oxide  volume  mixing  ratio  of 
-7 

5  x  10  in  the  upper  stratosphere  (Keneskea  1967).  The  model  atmosphere 
for  the  computations,  used  primarily  for  the  temperature  and  total  density 
profiles,  is  taken  from  the  U.  S.  Standard  Atmosphere  Supplements,  1966 
and  uses  the  30°  latitude  January  model  below  120  km  and  the  Spring/Fall 
model  with  exospheric  temperature  1000  K  above  120  km. 

The  results  shown  are  not  qualitatively  different  from  those  of 
Corbin,  et  al.  (1969).  The  quantitative  differences  result  primarily  from 
fee  use  of  different  molecular  abundances.  The  use  of  a  lower  exospheric 
temperature  has  negligible  effect  because  at  the  higher  altitudes  where  the 
temperature  is  lower,  the  radiation  is  determined  by  absorption  of  radiation 
from  the  lower  atmosphere  or  the  sun  and  temperature  dependent  collisional 
processes  play  a  negligible  role.  The  nitric  oxide  fundamental  band  at  5«  Zjjun* 
is  an  exception,  because  of  the  excitation  by  collisions  with  atomic  oxygen. 
Even  for  this  molecule,  results  are  not  appreciably  different  during  the  day. 
At  night,  however,  limb  radiances  at  5.  3  yxn  do  not  fall  off  as  rapidly  with 
altitude  below  200  km.  The  carbon  dioxide  bands  added  do  not  make  an 
appreciable  difference  in  the  limb  radiances  because  the  basic  excitation 
mechanism  for  these  bands  is  absorption  of  solar  or  terrestrial  radiation  in 
the  weak  combination  and  overtone  bands  at  5. 18  and  4.82yUm. 
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Limb  Viewing  Spectral  Radiances,  17  to  25 /nxn 
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The.three  general  areas  of  investigation  in  this  research  are  (a)  the 
chemistry  of  the  upper  atmosphere  and  its  relation  to  production  of  infrared 
radiation,  (b)  the  physical  mechanisms  wjhich  excite  and  quench  molecular 

I.  * 

quantum  levels  which  can  radiate  in  the  infrared,  ."vnd  (c)  improvement  of  models 

of  the  transport  of  infrared  radiation. 

*  1 .  » 

The  principal  results  of  the  current  program  have  been  the  development 
of  a  program  to  compute- the  abundances  of  infrared  emitting  molecular 

s 

species  insofar  as  knowledge  of  chemical  reaction  rates  permits,  and  an 
extension  of  the  limb  viewing  infrared  radiance  program  of  Corbin,  et  al.  (1969). 

,  This  work  is  incomplete,  but  permits  the  definition  of  areas  in  which  future 
work  will  be  profitable. 

i 

In  principle,  especially  at  the  higher  altitudes  considered  here,  the 
chemistry  and  diffusive  mechanisms  if  properly  described  should  permit  the 
predicted  concentrations  to  approach  a  diurnal  steady-state.  That  is,  the 
various  concentrations  should  be  reattained  in  every  twenty-four  hour  period. 
Likewise  there  should  be  other  steady-state  cycle  achieved,  namely,  those 
varying  with  the  seasons  and  the  11 -year  solar  cycle.  These  latter  effects 
require,  of  course,  much  more  elaborate  modeling  than  is  provided  in  the 
present  context  and  must  consider  variation  in  the  total  amount  of  solar 
energy  reaching  the  various  portions  of  the  earth's  atmosphere  and  the  large 
scale  ^ind  patterns  generated  by  these  variable  inputs. 

It  was  not  possible  in  the  present  program  to  exercise  it  sufficiently  to 
allow  the  diurnal  steady  state  behavior  to  be  investigated  over  the  desired 
range  of  altitudes,  species,  and  uncertain  reaction  rates.  One  run  through 
a  24-hour  cycle  was  performed  over  a  restricted  altitude  regime  to  determine 
how  much  the  solution  varipd  during  a  single  day  from  the  specified  initial 
conditions.  A  separate  run  was  pferformed  that  covered  a  total  of  12-days  real 
time.  These  results  were  not  acceptable  due  to  the  longer  time  steps  that 
were  required  in  the  program.  The  exercise  of  the  code  must  be  continued 
to  deterprine  the  permissible  time  steps  that  will  allow  steady-state  to  be 


'  reached  within  a  reasonable  machine  running  time.  Investigations  must  also 

be  made  of  the  effects  on  this  process  of  the  uncertainties  in  the  experimental 

j  reaction- rates. 

! 

j  At  altitudes  above  100  km,  the  principal  source  of  error  in  computed  limb 

* 

•  radiance  profiles  is  the  knowledge  of  molecular  abundances.  The  individual 

i 

|  lines  in  the  molecular  bands  are  weak  enough  that  the  neglect  of  collision 

i 

j  broadening  effects  which  lead  to  a  Lorentz  or  similar  shape  in  the  wings  of 

‘  lines  are  negligible.  At  lower  altitudes,  the  effects  of  the  line  wings  will 

j  increase  actual  radiances,  and  further  work  is  required  to  include  this  effect 

in  band  modeling.  The  effect  of  neglecting  the  far  wings  of  the  lines  is  expected 
to  be  greatest  in  the  carbon  dioxide  bands,  and  here  the  problem  is  complicated 
by  the  non-Loreutz  wing  shape.  Further  work  should  be  done  on  band  modeling. 

Measurements  of  radiation  in  spectral  regions  outside  the  nominal  limits 
five  to  twenty-five  microns  are  available  and  can  aid  in  defining  the  processes 
important  in  producing  radiation  in  that  spectral  region.  The  capabilities  of 
the  limb  viewing  atmospheric  radiance  model  should  be  expanded  to  include 
the  carbon  dioxide  and  water  vapor  bands  at  Z.7  micrometers  and  other  known 
bands  for  which  sufficient  data  exists  in  the  spectral  range  between  2.7  and  25 
micrometers.  Estimates  of  global  variations  of  infrared  radiation  are  required 
and  will  require  at  least  a  qualitative  estimate  of  the  effects  of  horizontal 
transport  (winds,  global  circulation  patterns)  on  chemistry. 

Further  work  suggested  may  be  summarized  in  the  context  of  the  present 
problem  areas: 

Chemistry.  The  main  effort  here  has  been  to  develop  an  efficient  and 
reliable  program  to  integrate  the  photochemical  reaction  rate  equations, 
including  vertical  transport  by  eddy  mixing  and  molecular  diffusion,  for 
atmospheric  atoms  and  molecules  which  radiate  in  the  infrared.  Lt  present 
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eighteen  different  neutral  atoms  and  molecules  are  included.  These  are 
molecular  nitrogen,  N£,  which  is  treated  as  a  background  gas  with 

concentrations  determined  only  by  tiansport;  two  states  of  atomic  nitrogen, 

4  2  3  1 

N(  S)  and  N(  D);  two  states  of  atomic  oxygen,  0(  P)  and  0(  D);  molecular 

oxygen,  O^;  ozone,  O^;  the  hydroxyl  radical,  OH;  the  hydroperoxyl  radical, 

HO  ;  water  vapor,  H  O;  hydrogen  peroxide,  H-O  ;  atomic  and  molecular 

C*  Ct  Ct 

hydrogen,  H  and  H^;  nitric  oxide,  NO;  nitrogen  dioxide,  NO^;  methane, 

CH^;  carbon  dioxide,  CO^;  and  carbon  monoxide,  CO.  In  addition,  argon 
is  included  to  serve  as  a  further  check  on  transport  calculations.  Running 
time  is  slower  than  desired,  and  further  development  should  include  a  search 
for  a  more  efficient  integration  scheme,  as  well  as  the  necessary  addition  of 
ionized  chemical  species.  Qualitative  estimates  of  the  effects  of  horizontal 
transport  on  atmospheric  chemistry  should  be  made. 

Excitation  Mechanisms.  Data  on  excitation  of  vibrational  levels 
of  infrared  emitting  molecules  by  collision  with  neutral  molecules,  ionized 
molecules  and  electrons  need  to  be  continually  reviewed.  Excitation  of 
vibrational  and  electronic  states  of  other  molecules,  particularly  nitrogen 
and  metastable  states  of  atomic  and  molecular  oxygen,  which  can  transfer 
energy  to  infrared  emitting  molecules  should  be  included.  As  relevant  new 
rates  are  found,  they  should  be  added  to  the  chemistry  and  radiance  computer 
programs. 

Radiative  Transport  Modeling.  The  addition  of  new  bands  requires 
further  modeling  of  individual  bands.  The  current  method  should  be  improved 
and  compared  with  other  methods. 
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RETURN  SINT  40 

0  CONTINUE  SINT  50 
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